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ABSTRACT 


The purpose of this paper is to define and present the wicking material properties 
that are considered to be important to the operation of a vapor-chamber fin or 
a heat pipe. 

Three classes of wicking materials are studied: sintered metal screens, sintered 
metal powders, and sintered metal fibers. The porosity of these sintered ma- 
terials ranges from 47. 7 to 91. 8 per cent. 

Two characteristics of the wicking material are considered to limit the operation 
of the heat pipe, 1) the capillary pumping characteristics of a wick, and 2) the 
evaporative heat transfer characteristics of a liquid-saturated wick. 

To evaluate the effect of capillary pumping characteristics, a simplified analysis 
of a planar-wick heat pipe is made. The result of this analysis gives the maximum 
operating length of a planar-wick model in terms of the external boundary con- 
ditions, the heat pipe fluid properties, and the capillary pump characteristics 
of the wicking material. The capillary pump characteristics are found to be 
proportional to the equilibrium height to which the heat pipe liquid will rise in 
the wicking material, divided by the wicking material friction factor, when net 
body forces are zero. The latter is the reciprocal of the permeability for the 
wicking material. When gravity effects cannot be neglected, then the equilibrium 
height to which the heat pipe liquid will rise in the wicking material becomes 
another distinct factor. 

The results of wick equilibrium height experiments and wick permeability experi- 
ments run on the three classes of wicking materials are presented. Both water 
and Freon 113 are used as working fluids in these experiments. These results 
are combined to yield the capillary pumping characteristics of wicking materials. 

To evaluate the effect of the evaporative heat transfer characteristics of wicks, 
experimental data on porous samples selected from the three classes of wicking 
materials is presented. This data results from evaporative heat transfer 
experiments run on planar-wick samples using water and Freon 113. Compari- 
son of results with the heat transfer characteristics of a flat plate submerged 
in water indicate that equivalent or superior performance can be obtained with 
wick covered surfaces. However, the data also indicates that the entrapment 
of vapor bubbles in the wick matrix may cause premature film boiling in the 
porous material at relatively low heat fluxes , depending on the structure of the 
wicking material. The heat transfer characteristics of Freon 113 were poorer 
than those exhibited by water. 
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I. SUMMARY 


This report summarizes the work performed to obtain the capillary pumping and 
evaporative heat transfer characteristics of wicks with respect to vapor-chamber 
fin or heat pipe operation. 

To determine the capillary pumping characteristics which affect the operation of a 
vapor-chamber fin, a simplified analysis of a planar-wick geometry was made. 

The simplified analysis indicates that the maximum height to which a liquid will 
rise in a vertical wick, j2 m , and the wick friction factor, Ki, are important para- 
meters in deter ming the capability of a capillary pump in a heat pipe. The larger 
the value of the capillary pump parameter, the better the capillary pump, 

when head effects due to gravity are negligible. If an application involves opera- 
tion at a condition in which a capillary pump must work against gravity, considera- 
tion must be given to both the value of Jt m /K^ and ofj£ m itself. The important fluid 
properties for capillary pumping are also indicated by the simplified analysis. 

The equilibrium height of liquid rise in a vertical wick was experimentally deter- 
mined for 10 different wick samples with water as the test fluid. Equilibrium 
heights of six samples were also determined with Freon 113. Three classes 
of wicking materials were studied: sintered metal screens, sintered metal 
powders, and sintered metal fibers. The porosity of these sintered materials 
ranged from 47. 7 to 91. 8 per cent. The wick friction factor was determined for 
thirteen different samples with water as the test fluid. Tests were repeated on 
four of these samples with Freon 113 as the test fluid. The friction factor of 
each sample was found to be essentially constant over a range of mass velocities 
and temperatures for the two fluids tested when they were deaerated. When the 
test fluid was aerated water, the wick friction factor was found to increase with 
time. 

The experimentally determined capillary pumping characteristics of the three 
classes of materials indicated that as a class the sintered fiber samples were 
best, followed by the sintered powder samples. The sintered screen samples 
demonstrated the lowest potential as capillary pumps. 

The evaporative heat transfer characteristics of six different samples from the 
three general classes were determined experimentally. The effects of orientation, 
thickness of sample, degree of submergence of the sample in the test fluid, test 
fluid and wick surface modifications were determined. Some of these porous 
materials provide high values of heat transfer coefficients and limiting heat trans- 
fer rates for boiling when adequate paths for vapor escape from within the wick 
are provided. Other samples which do not provide vapor removal paths result 
in values of heat transfer coefficients significantly lower than those of a bare 
flat surface. Three sintered fiber samples were determined to have the best 
heat transfer characteristics of those tested, followed by a sintered powder 
sample and sintered screens. 
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II. INTRODUCTION 


1 2 3 

The work of Grover, et al , and other investigators * has shown that the heat 
pipe is a heat transfer device that can exhibit an extremely high effective thermal 
conductivity, much greater in fact than any known homogeneous material. The 
heat pipe consists of a long closed container in which vaporization and condensa- 
tion of a fluid take place. Heat added at one end of the container causes evapora- 
tion of liquid into vapor. Condensation of the vapor along the length of the con- 
tainer maintains the surface at a nearly constant temperature. The resulting 
condensate is returned to the heated end of the container by the action of capillary 
forces in the liquid layer which is contained in a wick lining the inside of the 
cavity. 

4 

A parametric study done by Haller, Lieblein, and Lindow indicates that the heat 
pipe might be used as a "vapor-chamber fin" in reducing the weight of a radiator 
for a Rankine-cycle space poweiplant. Such a fin is shown in cross-section in 
Figure 1. Possible vapor-chamber fin configurations for use in a space radiator 
are shown in Reference 4. 

An investigation was begun under Contract NAS 3-7622 to explore and define the 
mechanisms of fluid and heat transport in vapor-chamber fins or heat pipes, to 
provide design information for space radiators and other applications. The 
investigation is divided into three parts, 1) wicking studies, 2) boiling studies, 
and 3) steady- state fin operation. This report summarizes the accomplishments 
of the first two parts of the program. 

In order to obtain an understanding of fin operation, wicking and boiling studies 
were conducted in order to evaluate the capillary and frictional forces associated 
with the wicking materials and the boiling heat transfer characteristics of wicks 
to be used in the fin. The mechanistic model and simplified analysis of fin opera- 
tion presented in Section HI of this report indicate the need for this information. 
Section IV lists and explains the selection of wicking materials and fin fluids which 
were tested. The wicking material studies are treated in Sections V, VI and VII, 
where an explanation and description of the equations and experimental apparatus 
necessary to evaluate the forces acting on a liquid flowing in a wick are given. 

The test program to measure the boiling and evaporation heat transfer properties 
of a liquid in a wicking material is discussed in Section VUI. 


See Appendix 1 for numbered references 
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r— LIQUID RETURNED TO EVAPORATOR BY CAPILLARY FORCES 


EVAPORATING 

LIQUID 


REJECTED HEAT 


CONDENSING VAPOR 


LAYER OF 

WICKING 

MATERIAL 


Figure 1 Vapor- Chamber Fin in Operation 



m. ANALYSIS OF VAPOR-CHAMBER FIN 


A. Vapor-Chamber Fin Operation 

The various mechanisms that control the operation of the vapor-chamber fin 
can be discussed by considering the schematic diagram of the fin shown in 
Figure 2. This diagram illustrates the fact that the vapor-chamber fin is similar 
to a heat engine consisting of an evaporator, a condenser, and a pump. Although 
there is no net work output, pumping work is provided by capillary forces in 
transporting the liquid from the condenser to the evaporator. 

In the evaporator section of the fin, heat must be conducted from the source of 
heat through the metal wall of the fin and the composite layer of liquid and wick- 
ing material in contact with the wall. This energy transfer results in evapora- 
tion of the liquid at the liquid-vapor interface. 


LIQUID FLOW 



PAGE NO. 


4 



In the condenser section of the fin the reverse process to that in the evaporator 
occurs. Vapor condenses at the liquid-vapor interface, liberating its latent heat 
of vaporization. This energy is then conducted through the composite liquid-wick 
layer to the fin wall, then through the fin wall itself, from which it is rejected to 
the environment of the fin. Radiation is the mode of heat rejection from the fin 
in a space-radiator application. 

The capillary pump serves the function of returning the condensed liquid to the 
evaporator. The pump may consist of a wicking material, the pores of which 
are small enough to provide capillary forces of sufficient magnitude to cause 
the liquid to flow back to the evaporator. The differences in the radius of curva- 
ture of the liquid-vapor interface in the condensing and evaporating sections of 
the fin provide the driving force of the pump. The details of this capillary pump- 
ing will be discussed in Section m C below. 

When the various transport processes that occur during fin operation are analyzed, 
it is found that there are essentially two possibly interrelated phenomena which 
limit the maximum rate of energy that the fin can transfer. They are: 

1) The heat flux at which apparent film-boiling occurs in the wick of the 
evaporator, and 

2) The maximum liquid flow rate that the capillary pump can provide to 
the evaporator. 

There may be other limiting phenomena, depending on the particular application 
in which the fin is used. For example, the thermal conductivity of the liquid-wick 
layer may be of such a magnitude as to provide a significant resistance to heat 
flow. But for the present study in which conventional liquids and space-radiator 
geometries will be studied, the film-boiling heat flux and the maximum capillary 
pump output seem to be the two limitations of most importance. Thus, although 
other factors that influence the operation of the fin will be taken into account in 
this program, these two phenomena will be studied separately and in greater 
detail. 

It is unfortunate that these two limiting phenomena depend on the wick; for in 
the near zero-gravity environment of a space radiator, it is conceivable that 
a vapor-chamber fin wick would not be needed. The liquid would simply wet 
the walls of the fin cavity and due to changes in the radius of curvature at the 
liquid-vapor interface, would flow back to the evaporator section of the fin. In 
the present study all tests will be run in a 1-g environment, thereby necessitat- 
ing the use of a wicking material in any fin models tested in order to make the 
capillary forces large relative to the gravitational force. However, even in an 
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actual vapor-chamber fin space radiator, a wick of some form would probably 
be used to ensure that the radiator would work in a force field such as gravity. 

In all discussions that follow it is assumed that the fin is in a gravitational field. 

When either of the two limiting phenomena occurs, the result should be the same: 
any attempt to increase the net heat input to the fin would result in a large tem- 
perature rise at the base of the fin. This temperature increase, shown in the plot 
of Figure 3, is probably caused by an insufficiency of liquid at the liquid-wick- 
solid interface in the evaporator of the fin. 



Figure 3 Idealized Performance of Vapor-Chamber Fin as Function of Heat 
Input 

Thus, when the performance of a vapor-chamber fin falls off as the heat input 
is increased, the question arises as to which of the two limiting phenomena is 
responsible, wick film-boiling or insufficient capillary pump output. In order 
to answer the question a separate evaluation of each of these two possible causes 
of fin failure is needed. The following two sections discuss each possibility. 
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B. Limitation Due to Wick Boiling 


An evaluation of the first of the limiting phenomena listed in the preceding section 
involves a study of the process of boiling of a liquid in a wick. 

The actual process occurring in the evaporator section of the fin can be visualized 
in the following manner. As heat is initially added in the evaporation region, the 
heat will be conducted through the fin wall and the wick-liquid composite to the 
liquid-vapor interface where evaporation will take place. As the heat input to 
an operating fin is increased, it is presumed that the liquid- wick- solid interface 
temperature will rise to the point where the adjacent liquid is sufficiently super- 
heated to cause nucleate boiling at the interface. If an equilibrium condition is 
to occur without excessive superheat, the capillary and/or buoyancy forces must 
be great enough to cause convection of the resulting vapor bubbles through the 
wick to the vapor cavity of the fin and supply fresh liquid to the heating surface 
at a sufficient rate. The buoyancy forces will be dependent upon the gravitational 
field in which the fin is operating. If vapor does not leave the interface at the 
rate at which it is produced, liquid will be excluded from the interface, causing 
film-boiling to set in. The effects of such a process are shown in Figure 3, as 
a sharp increase in temperature at the base of the fin, causing the fin to eventually 
become ineffective. 

Because of the limited success with theoretical models of the simpler case of 
pool boiling, an attempt to evaluate boiling heat flux limits analytically in a 
relatively complicated wick structure appears unpromising. Thus, an experi- 
mental program to study boiling in wicks was conducted and will be described 
in Section VIII. 

C. Limitation Due to Capillary Pump 

The second major limitation on fin operation is the maximum output of the 
capillary pump. Stated differently, it is a matter of the wicking forces of the 
wick material overcoming the frictional and gravitational forces to provide 
sufficient flow rate. 

To analyze the action of the wick as a capillary pump, the simplified model of 
the vapor-chamber fin shown in Figure 4 can be considered. The liquid pres- 
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LIQUID-VAPOR INTERFACE 


VAPOR CAVITY 



Figure 4 Two-Dimensional Model of Vapor-Chamber Fin 


sures at the ends of the fin are (Pl)^ and (Pl) 2> while the vapor pressures at 
the ends directly above the liquid are (Py)^ and (Py )2 respectively. In all of 
the following discussions it is assumed that no boiling takes place. 

For the fin to operate, there must be a liquid flow from 2 to 1. Therefore, to 
overcome frictional forces and momentum changes 


< P L>2 


< P L>1 


Similarly, there must be a counter vapor flow from 1 to 2, so that 


( 1 ) 


( p v >! S < P V >2 < 2 > 

If the system in Figure 4 is at steady state, the sum of the pressure changes 
that a given fluid particle experiences as it leaves the evaporator and returns 
again, is zero, or 

[< p v)i-< p v> 2 ] + [<V2-< P L> 2 ] + [< P L>2 - < P L>1 ] + 


[< p L >r<Vi] =° 


( 3 ) 
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As the liquid flows in the wick or the gas in the cavity, the pressure changes are 
balanced by frictional and inertial forces. These pressure drops due to friction, 
rate of change of momentum, etc. , can be defined as 

(AP) d s ( Ap y) d + ( AP L) d s <Vl ~ + ^ 1?2 ~ ^ P L^1 ^ 

where, in general: 


(AP y ) d > 0 and (AP^ > 0 


( 5 ) 


These terms are usually positive since they consist primarily of the frictional 
components. 


At any vapor-liquid interface there is some radius of curvature of the resulting 
surface, whether it be a very small value or infinity. The Laplace-Young equa- 
tion of capillarity 5 shows that there is a pressure change across an interface, 
given by 


AP = 



( 6 ) 


where AP is the pressure difference across the surface, o is the liquid-vapor 
surface tension, and R and R' are the two radii of curvature necessary to describe 
any three-dimensional surface. The highest pressure is always on the concave 
side of the interface when R and R T have identical signs. When it is not clear 
which is the concave side (e.g. , a saddle-shaped interface) the higher pressure 
is on the side of the smaller radius. 


If the liquid wets the wick, that is, if the contact angle is less than 90°, and if it 
is assumed that the condenser and evaporator vapor-liquid interfaces can each be 
characterized by a single radius of curvature, the remaining terms in Equation 
(3) become 


< P V>2 - < P L>2 


2 a 


R„ 


( 7 ) 


< P L>1 - < P V>1 


- 2a 

R- 
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where <r is the liquid-vapor surface tension and and R 2 are the radii of 
curvature describing the liquid- vapor interfaces in the evaporator and condenser 
regions. 

Since the fin operates at a nearly constant temperature, any temperature -dependent 
variation of * can be neglected. Combining Equations (3), (4), (5) and (7) yields 

2 <r 2 <r 

< AP >.r -if- - — > 0 < 8 > 

This means that 


Thus, the capillary pump, that is, the wick, overcomes the frictional and inertial 
forces in the vapor-chamber fin by allowing a larger radius of curvature to form 
in the condenser than in the evaporator. Of course, the liquid must wet the wick 
for this to be true. 

Physically this means that vapor condensing out on the wick causes a large inter- 
facial radius of curvature to form as the liquid tends to flood the wick. However, 
the heat added to the evaporator makes the vapor-liquid interface retreat into the 
matrix of the wick, causing a smaller interfacial radius of curvature, a radius 
probably characteristic of the size of the pores in the wick itself and the contact 
angle. These differences in radii of curvature at different locations along the 
wick make the porous medium of the wick act as a capillary pump. 

At first it might be thought that to increase the capacity of the pump it would be 
best to decrease the size of the pores in the wick, thereby decreasing R^ in 
Equation (8). However, this also tends to increase the frictional forces in the 
wick. The relationship between the capillary forces, frictional forces, and in- 
ertial forces will be derived in the next section. 

The upper limit on the frictional and inertial pressure drops in both the vapor 
and liquid can be calculated using the above expressions. The upper limit on the 
pressure drop through the vapor is obtained by using Equations (4), (5) and (8) 
and assuming that the liquid pressure drop is zero to get the inequality 

<Vi-< p v>2 < 2<r (~iq- - “ir) (9) 
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The right side is maximized when the condenser is completely flooded, that is, 
when R 2 — > cd. Hence, for the fin to operate 


< P V>1 - < P V> 2 < < 10) 

Values of the right-hand side can be obtained by considering the minimum 
effective pore size in the wick, which for most wicking materials considered 
yields 


— m 5 to 20 inches of H.O 

2 


By similar reasoning it can be shown that the maximum pressure drop in the 
liquid is bounded by this same value. 

It has been demonstrated that the shape of the liquid-vapor interface is critical 
to the operation of the capillary pump. In particular, the radii of curvature in 
the evaporator section are most crucial, for as long as the evaporator supplies 
sufficient vapor and the fin is sufficiently charged with liquid, the radii of curva- 
ture in the condenser will be large, making the pump term in Equation (8) strongly 
dependent on R^. 

It is not apparent what the actual shape or position of the entire liquid-vapor 
interface is in the fin evaporator. Three possible positions are shown in Figure 
5. 


The situation pictured in (a), where the interface has retreated only slight- 
ly into the wick matrix, is a correct representation for lower heat-input 
rates (assuming that the wick has been initially saturated with liquid). However, 
the situation cannot be so easily visualized while the wick is providing liquid flows 
at higher heat-input rates, for under such conditions, the liquid-vapor interface 
is certain to retreat further into the wick. At still higher values of heat input, 
the wick will no longer be able to supply sufficient liquid to be evaporated. The 
wick matrix will then dry and the interface could assume the position shown in 
(b) or (c), or be a combination of these. Whatever the interfacial position, when 
the maximum flow rate of the pump is exceeded the evaporator will dry and the 
base temperature will rise, as shown in Figure 3, causing eventual loss of fin 
effectiveness. 
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Because it is not clear what the shape and position of the evaporator liquid-vapor 
interface will be at high heat-input rates, it is essential that an isotropic wicking 
material be used in the present study of the vapor-chamber fin (the term isotropic 
will be clarified in Section IV A). In this way any limitation in fin operation due 
to the maximum flow rate provided by the capillary pump cannot be attributed to 
the directional properties of the wick. 



LIQUID-VAPOR 

INTERFACE 


LIQUID-VAPOR 

INTERFACE 


LIQUID -VAPOR 
INTERFACE 


Figure 5 Simplified Picture of the Liquid-Vapor Interface in the Evaporator 
Section of the Fin 

So far it has been shown how the capillary pump functions qualitatively. In the 
next section a simplified analytical model of the fin will be derived, based 
entirely on the limitations imposed by the capillary pump. For the present it 
will be assumed that the evaporator section is sufficiently small compared to 
the condenser section so that the pressure drop of the liquid in the evaporator 
section may be neglected. It will also be assumed that the evaporator heat flux 
is below the film boiling value. 

With the help of some simplifying assumptions, an expression will be derived 
that will show what effect the capillary pump capability has on the maximum 
size of a vapor-chamber fin for a given heat input, or conversely, what effect 
the pump has on the maximum heat input that a fin of a given size can accept. 
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D. Simplified Model of Vapor-Chamber Fin 


This section contains the derivation of a simplified model of the vapor-chamber 
fin. An analytical model of the heat pipe has been formulated by Cotter 6 . The 
purpose of this model of the vapor-chamber fin is to discover the parameters 
necessary to adequately describe the mechanism of the fin. Within the frame- 
work of the assumptions made in its derivation, the model includes the effects 
of rate of change of momentum, surface tension forces, and the frictional forces 
in the body of the wicking material. 



Figure 6 Simplified Model of Vapor- Chamber Fin 


Figure 6 shows a simplified model of the vapor-chamber fin that is rejecting 
heat at a specified heat flux density q. The fin width in the z-direction is b and 
the heat rejecting length is x m . Because limiting boiling heat fluxes are not 
being considered, the evaporator section, from x = x m to x = xrp, is assumed 
to be of sufficient size to provide enough vapor to meet the requirements set by 
the specified value of the condenser rejection flux density q. In addition, the 
following assumptions are made concerning the model shown in Figure 6 ; 

1) The wicking material is rigid, is of constant thickness 8 , and is 
saturated with a wetting liquid along the entire condenser length x m . 
Also the wicking material is considered to be isotropic, so that at 
any given cross-section of the wick, the pore area, A p , and the total 
area of the cross-section, A rp, are related by 

A 

-r 2 - = • (ii) 
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where « is the porosity of the material. This relation is derived in 
Section IV B. 

2) The vapor and the liquid are at the same constant temperature from 
x = 0 to x = x m . This assumes no subcooling of the liquid. It 
follows that the vapor pressure Py is also constant along this length, 
and that there is little conducted energy or variation in convected 
energy in the x or z directions within the wicking material. 

3) The vapor that condenses on the liquid-vapor interface has a velocity 
Uy in the y-direction. This velocity has no x or z components; thus 
there are no momentum contributions in these directions. 

4) The effects of gravity are not considered. 

5) At any given point, the liquid-vapor interface can be characterized by 
one radius of curvature, R. 

6) The velocity of the liquid in the wicking material, u^, has an x com- 
ponent only, and is uniform across the wick thickness 5 . This 
velocity is considered to be an average velocity in the pores of the 
wick. 

7) Any terms containing differentials raised to the second or higher powers 
are neglected. 

8) The heat flux is considered constant with length. 

Figure 7 shows the forces and momentum fluxes associated with a differential 
element of the saturated wicking material at a given x in the heat rejecting 
portion of the fin. 


Applying the continuity equation for the element of Figure 7, the relationship 
between the liquid flow rate w^ and the vapor flow rate Wy is obtained 


Because 


and 


(W L>x + W V = ( Vx + dx 


(12) 

(W L ) x = P L ' < 5 b > U L 


(13) 

/ 

dU T \ 


< W L>x + dx =P L‘ <« b >K + 

ax L H 

(14) 
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PRESSURE = P v 



Figure 7 Differential Element of Wicking Material for Momentum Equation 


where p ^ is the liquid density, it follows that 



dU L 



W V = P L dx d **<* b > 

(15) 

where 

w v = p v “v (bdx) 

(16) 


and p y is the vapor density. 
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The momentum equation in the x-direction for the element in Figure 7 is 


2F = M , 
x x + dx 


- M 


x 


( 17 ) 


The momentum terms are 

Pj 2 

M = — — u^ «(b 5) (18) 

x g Q L 

2 

P T 9 Pj d(u ) 

M = — u! e (b5) + — — ■ dx c (b $ ) (19) 

x + dx g L g dx 


where g is the proportionality constant in Newton’s second law. It should be 
noted that because of Assumption 3 there is no momentum contribution to 
Equation (17) by the vapor as it enters the top of the element. Calculations show 
that when such a momentum term is included, it is negligible compared to the 
other terms in Equation (17). In addition, the presence of the rigid wick at the 
liquid- vapor interface interferes with the interaction of the liquid and vapor 
with respect to momentum, and thus makes the inclusion of such a term doubtful 
in the analysis of an actual vapor-chamber fin. 


Two forces make up the force term in Equation (17): the pressure force due to 
the differences in radius of curvature at x and x + dx, and the frictional force 
Ff acting on the liquid as it flows through the porous medium of the wick. 

Thus 


SF = ( F 
x p 


>x-< 


F ) - F 

p'x + dx f 


( 20 ) 


The pressure force Fp is obtained by use of the Laplace- Young equation, Equa- 
tion (6). Because of Assumption 5 this reduces to 

2 <r 

AP = -g— (21) 

where AP is the pressure drop across the liquid-vapor interface at x, a is the 
liquid-vapor surface tension, and R is the radius of curvature of the interface 
at x. Thus the pressure forces at x and x + dx are 


< F p>x = ( P V * <* b) 


( 22 ) 
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< F P>x + dx= ( F v- RT5T ) « < 8b> < 23 > 

where the vapor pressure Py is constant, by Assumption 2. 

The net pressure force on the element is found by adding Equations (22) and 
(23) vectorially to get 


d"R 

(F P>x-< F P>x + dx = - ( < 5b »< 2 ') — r < 24 > 

R 

where RdR has been neglected when it is compared to R^. 

Because of the low flow rates and velocities encountered with capillary flow in 
wicks, the flow produced by the net capillary force in Equation (24) can be 
assumed to be both laminar and relatively free from inertia effects. As pointed 
out by Scheidegger 7 , both of these conditions must be met if Darcy's law-' 7 for 
flow through porous media is to apply. Darcy's law can be stated as 



(25) 


where 



w L 


= pressure gradient in the porous medium due to 
frictional forces 
= liquid viscosity 
= liquid density 

?= total cross-sectional area of flow (i. e. , the sum of 
the pore area and the area of the solid at a given cross- 
section through the porous material) 

= liquid mass flow rate 


The term is a constant of proportionality. This constant is particularly 
important in the evaluation of wicking materials. It is dependent on the geometry 
and dimensions of the passages in the wicking material and could possibly be 
described adequately if an equivalent passage diameter, the wick tortuosity, and 
porosity of the material were all known functions. Since this is usually not the 
case, can be determined experimentally by passing a known flow rate of liquid 
through the wicking material and measuring the resulting pressure drop across 
the sample. Then Kp which will be called the wick friction factor, can be cal- 
culated from Equation (25). The constant has the dimensions of (length) -2 
and is the reciprocal of permeability. 
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The frictional force on the fluid element shown in Figure 7 is given by 


F f = ("fr) ^ 4 b >i < 26 > 

By combining Equations (13), (25), and (26), an expression for the friction term 
in Equation (20) is given by 

2 

F = K € ( 5 b) p L u L dx (27) 


Combining Equations (17), (18), (19), (20) (24), 
momentum equation 

dR 

- 2 ' “1“ ' * K i "l u l * = 


and (27) results in the differential 



d (u L ) 2 

dx 


dx 


(28) 


The remaining relationship necessary to describe the operation of the fin is the 
energy equation. Figure 8 shows the energy fluxes associated with a differential 
wick element at some x in the condenser portion of the fin, where Q is the energy 



Figure 8 Differential Element of Wicking Material for Energy Equation 


PAGE NO. 18 



transferred per unit time, the superscripts c and k refer to convection and con- 
duction, and L and V refer to liquid and vapor. 


Neglecting changes in kinetic energy, the individual terms in Figure 8 are given 
by 

out 


in 


c lc 
Q = Q + Q 


«L>x = W L h L 


<Q l ) = -k L ( 8 b) 


< Q L>x + ' = w ~ K + 


dx L L 


(Q.) . =-k T (6 b) 

^I/x + dx L 


dT. 


dx 


d < w L h L ) 

dx 

dT T 


dx 


[ 


dx 


d 2 T. 


dx 


dx 


] 


(29) 


Q = q (b dx) 


where h is enthalpy and k^ is the thermal conductivity of the liquid-wick com- 
bination. 

Because of Assumption 2 all of the conduction terms in Equation (29) are zero. 
Thus the energy equation for the element in Figure 8 is 


Q V + < Q L>x - < Q l>x + dx + Q 


(30) 


The energy entering the top surface of the element is 

% = <4 = w v\ (31) 

Combining this with the continuity equation. Equation (15), results in 

du 

= 1^ p L dx e ( 5 b) (32) 


Putting Equations (32), (29) and (13) into Equation (30) yields the combined 
energy -continuity differential equation 
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dh. 


U L P L ' <\ - h L> P L : 


* i 
dx 


5 € 


= 0 


(33) 


But for an incompressible liquid where internal energy is a function of temperature 
only, the enthalpy variation with x is expressed as 


dx 


dT 


dP, 


= C 


L dx 


dx 


(34) 


where is the specific heat of the liquid. As shown in Section III C , the 
pressure drop in the liquid is very small, so that 


dP, 


dx 


~ 0 


(35) 


Also by Assumption 2 


dT 


dx 


= 0 


(36) 


Thus 


dx 


(37) 


so that Equation (33) becomes 

du 


* V * • P L 


(38) 


where h^ L is the latent heat of vaporization of the liquid. 

Integrating Equation (38) and using Assumption 8 yields the combined energy- 
continuity equation 


u = 


L *Vl { f P L 


x 


(39) 
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I 


showing that the liquid velocity varies linearly with x. The liquid flow rate then 
becomes 


w 


L 


__qb 

V 


x 


(40) 


Finally, Equations (39) and (28) can be combined to give an energy-continuity- 
momentum integral equation 





x dx = 


/. 


m 


2q x dx 


o*Vl p : 


(41) 


where, as shown in Figure 6, x m is the length of the heat rejection portion of 
the fin. 


To evaluate the maximum output of the capillary pump, Equation (40) shows that 
x m must be a maximum. In Equation (41) it is seen that to maximize x m the 
limits of integration on the surface tension must contain both minimum and maxi- 
mum values of R. In Section HI C it was shown that at x = 0, R can be very 
large and can be considered to approach an infinite value for optimum fin opera- 
tion. Similarly, as shown by Equation (8) the radius of curvature in the evaporator 
must be at some minimum radius of curvature, R min . To calculate the maximum 
value of x m , the assumption is made that at 

x = 0, R -» co 


and that at 


x = x 


m’ 


R 


R min 


(42) 


where R m i n is the minimum radius of curvature that the wicking material can 
support. It is further assumed that the entire evaporator liquid-vapor interface 
can be characterized by although no assumption has to be made about the 

position of the interface in the evaporator. 


Using the above limits for R, Equation (41) can be integrated and solved for x m 
to obtain 

72 


x - 
m 


P L 




R 


mm 


qKj_ / 

2 « \ 


1 + 




2a V 




(43) 
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If a range of representative values is substituted for the quantities in the non-unity 
term in the denominator of Equation (43), it is found that 




g o P L Vl 5 e2K i 


« 1 


(44) 


This term represents the contribution of the rate of change of liquid momentum, 
and can be neglected when compared to unity. Thus the expression for the maxi- 
mum length of the heat rejection portion of the fin is given by 


x 

m 



(45) 


In Equation (45) the liquid properties enter as the dimensional parameter 


N 


p l Nri/ 


(46) 


To maximize the performance of the fin, a liquid should have as large a value 
of N as possible. N for conventional liquids ranges roughly from 1 to 350 
Btu hr lb f /ft 3 lb m . By way of contrast, the liquid metal bismuth has a value 
of N = 2097 Btu hr Ibf/ft^ lb m at 1300°F. Other values of N are listed in 
Section IV D. 


The only quantity in Equation (45) that is not readily measurable is R m i n * Be- 
cause a maximum value Of x m is desired, R m in * s assumed to be the smallest 
radius of curvature that the wick can support. The method used to determine 
an effective minimum radius of curvature in a wicking material was to observe 
the maximum height to which the liquid rises in a vertical sample of the wicking 
material. If the pores in the vertical wick through which the liquid rises can 
be thought of as single capillary tubes of an equivalent diameter, the effective 
minimum radius of curvature is 

2 g o* 

R . = r- (47) 

min 

l_i m 

where is the maximum height of the liquid in the vertical wick when exposed 
to a gravity field, g. Equation (47) is obtained from a static force balance on 
the liquid in a capillary tube. See Section V A for the derivation of Equation 
(47). R m in in Equation (47) is termed the effective minimum radius of curvature 
because in actuality the smaller pores in a wicking material are not single 
capillaries but are interconnected with pores of larger size. However, 
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Equation (47) should still be an effective measure of for a vapor-chamber 

fin application. This will be demonstrated in Section V A. 

By combining Equations (47) and (45) it is seen that the final expression for the 
maximum length of the condenser section of the fin is given by 


Vi 


O / P L Nl' N 

f—)t 

8 P L 

B 

[ \ ^ J 

[ q )' 

e o ° 

' WR K 1 J 


(48) 


where 


<r = surface tension 
g = acceleration of gravity 
g 0 = Newton’s second law proportionality constant 
PL = liquid density 

hyL = liquid latent heat of vaporization 
8 = liquid - wick thickness (constant) 

An = equilibrium height of liquid rise in vertical wick 
Ml = liquid viscosity 
q = specified rejection heat flux 
= wick friction factor 


and where the subscript WR refers to wicking-rise test conditions. 

Thus, x m is maximized when the wicking parameter is large. It should 

be emphasized that gravity effects have not been considered in Equation (48) and 
that g, the acceleration of gravity, enters only through the separate process of 
determining R m i n from Equation (47). 


Equation (48) can be expressed in terms of the maximum total fin heat input 
Qmax. Thus, instead of solving for x m and specifying a heat rejection flux q, 
the problem can be restated to get 


Q 


max 


= 2 




(49) 


where b is the fin width in the z -direction. This is the maximum total heat input 

that can be transferred to a vapor-chamber fin with a condenser of length x m . 

A plot of Equation (48) is shown in Figure 9 for water at 350°F and a wick 
thickness of 8 = 0. 1 inch. This plot shows the maximum size of the heat re- 
jection portion of the fin of Figure 6 for various heat rejection rates 
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CONDENSER LENGTH , x m~ INCHES 



Figure 9 Condenser Length as Function of Condenser Heat Flux and 
Wick Properties Computed from Equation (48) 



as a function of the capillary pumping parameter £ m /K-p Again it must 
be remembered that neither gravity nor the evaporator section size are 
considered. 

If the effects of gravity are included in the analysis, then the maximum total heat 
input that can be transferred is given by the following equation (see Appendix 3 for 
the derivation). 


Q 

max 


/ 2p L Vlf ' 

\/ b 8 

-V-S- Pl \ 


)\ X 

/ \g <T ) 

\ > 
L ' L 

x m 


2p L 2 *Vl 

S /_ 

6b \ . 


g o \ 

l sin o 
/ 


WR 

1 In 

J K 1 


(50) 


In most applications of the vapor-chamber fin, the evaporator section is small 
in area compared to the condenser section and may be assumed to have no effect 
on the condenser, as long as the limiting heat flux is not reached. However, if it 
is inherent in a heat- rejection problem to have low evaporator heat flux, the 
evaporator may be of significant length when compared to the condenser. The 
following equation includes the effects of liquid pressure drop in the evaporator 
(see Appendix 3 for the derivation) and gravity. 



(51) 


The simplified analytic model of the vapor-chamber fin presented in this section 
illustrates the need for a separate wicking study to determine the value of ^ m /K f 
for various wicking materials. The results of such a study are outlined in the 
next sections. 
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E. Summary 


In the present study of the vapor-chamber fin for use in space-radiator applica- 
tions, the two key factors in the operation of the fin as a heat rejection device 
seem to be: 


1) The heat flux at which apparent film-boiling occurs in the wick of the 
evaporator section of the fin, and 

2) The maximum liquid flow rate that the capillary pump can provide to the 
evaporator section of the fin. 


Due to the complicated nature of the phenomenon of boiling itself, a separate 
experimental program was conducted to evaluate the limits imposed on fin opera- 
tion by the onset of boiling in wicks. 


To evaluate the effect of the capillary pump on fin operation, a simplified analytic 
model has been derived in which it was found that the output of the capillary pump 
depends on the liquid parameter 


N - 


P L *Vl 47 


(52) 


and the capillary pumping parameter 4^/K^ when net body forces are zero. 

Both of these parameters should be large to maximize the performance of the 
vapor-chamber fin. 

The capillary pumping parameter is the ratio of the maximum capillary forces to 
the frictional forces on a liquid flowing in a wick. The constant i m is a measure 
of the magnitude of the capillary forces associated with the liquid -wick combination. 
Physically, 2 m is the equilibrium height to which the liquid will rise in a vertical 
sample of the wick. 

The constant is a measure of the magnitude of the frictional forces associated 
with a flow of the liquid through the wick under conditions where Darcy f s law 
applies. This constant is obtained by measuring the pressure drop across a 
wicking material sample through which a given amount of liquid is flowing. 

Thus to study the effect of the capillary pump on the operation of the fin, a 
separate wicking materials study was conducted to measure values of and K-^. 
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rv. WICKING MATERIALS AND TEST FLUIDS 


A. Description of Wicking Materials 

A wick can be defined as a porous material whose pores are of such size that a 
liquid can be moved through it by the action of capillary forces. The pores must 
also be highly interconnected rather than separate and discrete. The pore sizes 
should fall into a range where a favorable balance exists between the capillary 
forces and the frictional forces associated with the motion of a liquid. 

The wicking materials considered in this study for use in vapor-chamber fins 
were narrowed down to three classes of metallic porous materials. They are: 

1) sintered screen materials 

2) sintered powder materials 

3) sintered fiber materials 

The reasons for selecting these materials, which are outlined below, include 
ease of fabrication, dimensional stability, isotropy and thermal conductivity. 

1) Wetting Properties 

For a porous material to act as a wick, it must be wetted by the liquid. 

With average care in eliminating contaminants , nickel forms a wetting 
surface with the conventional liquids used in this study. Nickel oxide is 
even more likely to form a smaller contact angle than pure nickel, be- 
cause of its higher surface energy (see Figure 38). 

Because stainless steel is also wetted by most liquids, it was employed 
as a second choice. However much greater care must be taken with a 
stainless steel surface to insure that no contaminants are present to 
prevent wetting. The need to keep stainless steel especially clean has 

Q 

been pointed out by Fox, Hare, and Zisman. 

2) Chemical Activity 

Nickel, stainless steel, and their oxides are inert with respect to the work- 
ing fluids of this study. 

3) Fin Fabrication 

Any sample tested may eventuaUy be installed in a fin model. All of the 
porous materials listed are capable of being attached to the metal surface 
of the fin cavity and are capable of being fabricated in a variety of shapes. 
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4) Dimensional Stability 


As shown by the analytic model of the last section, the wick dimensions 
must be known with some accuracy. All of the above materials are sufficiently 
rigid so that they will maintain their dimensions when saturated with a liquid. 
This is not true of soft fibrous and woven materials such as Fiberglas, textile 
fabrics, etc. The latter materials may either swell upon absorbing liquid or 
become compressed by capillary forces at a liquid- solid -vapor interface. 

Sample rigidity also allows the pore size distribution to be determined, since 
there are no experimental techniques for measuring this quantity in nonrigid 
porous materials. 

5) Isotropy 

As pointed out in Section in C , it is essential that an isotropic wicking 
material be used. By isotropic, it is meant that the wicking material has 
the same pore size distribution in any given direction. Thus, by using 
isotropic wicks, the operation of the capillary pump can be studied without 
complications caused by wick directional properties. 

The porous materials listed are essentially isotropic, with the obvious 
exception of the metallic screen layers. An analysis of the latter shows 
that for a square weave structure, the departure from isotropy is not great. 

6) Wick Thermal Conductivity 

The thermal conductivity in the liquid-wick layer is a factor in fin perform- 
ance. The effective thermal resistance of this layer is reduced if a wicking 
material of relatively high thermal conductivity is used. 

This is shown by considering the two limiting cases of heat conduction in a 
porous medium as presented by Gorring and Churchill 9. The lower limit on 
the effective conductivity of a two-phase porous medium is found by assuming 
that the two phases conduct heat in series. The resulting expression is 
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where k is the effective thermal conductivity, « is the porosity, and k s and 
k-j^ are solid and liquid conductivities respectively. To maximinize k, k s 
should be as large as possible. 

The upper limit on the effective heat conductivity would occur if both phases 
conducted heat in parallel. This yields 

(i) ■*-•(>- 

max 

Again, k s should be as large as possible to maximize k, the effective con- 
ductivity. 

Hence, from the standpoint of wick thermal conductivity, metal porous 
materials will be superior to other more conventional wicking materials, 
due to their inherently higher conductivities. 

The following is a detailed description of each class of porous material. 

1. Sintered Screen Materials 



The four sintered nickel screen samples that have been studied are listed in 
Table 1. 


TABLE 1 


Sintered Screen Samples 


Sample No. 

Mesh Size, 
(wires per inch) 

Wire Diameter, 
inches 

M7 

50 

0. 0090 

M8 

100 

0. 0045 

M9 

150 

0. 0030 

M10 

200 

0. 0022 


The sintered screen samples consist of layers of nickel screen that have been 
pressed against a nickel foil backing and sintered together at points of contact 
in a furnace with an ammonia atmosphere. The mesh referred to is the number 
of wires per linear inch in either the warp wire direction or the woof wire 
direction. The wire diameter given is that of both the warp and woof wires. 
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A typical screen sample (in this case 100 mesh) is shown in cross-section 
in Figure 10 and in plan in Figure 11. The solid structure at the bottom of 
the sample in Figure 10 is the nickel backup foil to which the screens are 
sintered. The porous material of this sample is 0. 0675 inch thick and the 
backup foil is 0. 0248 inch thick. These photomicrographs exhibit a number 
of features of sintered screens as porous materials. For instance, as the 
screens are stacked in Figure 10, it would probably be incorrect to charac- 
terize the pore size of the material based solely on the spacing of the wires 
in one layer of screen. The photomicrographs show that screen layers (at 
least in the two top layers) are stacked in such a way that alternate layers 
are displaced from one another by one-half a wire spacing. This means that 
wherever there is a square opening (or pore) formed by two warp and two 
woof wires in one screen layer, there is a crimped wire elbow protruding 
into the center of the opening from the screen layer immediately below the 
first. This results in a smaller pore size and a denser structure. Also 
noteworthy is the lack of any relatively large pores that are continuous 
from the backup foil to the free surface of the porous material. 

Comparison of Figures 10 and 11 also shows that sintered screen samples 
are not strictly isotropic. However, because of the interlocking of the 
screen layers, (i. e. , the displaced stacking discussed above) the departure 
from isotropy is not great. 

2. Sintered Powder Materials 

The six sintered nickel powder samples investigated are listed in Table 2. 

TABLE 2 

Sintered Powder Samples 

Powder Size Range, 


Sample No. 

Microns 

Ml 

150-297 

M2 

150-297 

M3 

297-841 

M4 

297-841 

M5 

297-841 

M6 

297-841 
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Mag:21X 

Figure 10 Cross-Sectional View of 100-Mesh Sintered Nickel Screen 
Sample M8 



rj Figure 11 Top View of 100-Mesh Sintered Nickel Screen Sample M8 
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These samples were made by compressing the nickel powder against a nickel 
backup foil, using a die. This "green” material was then placed in a furnace 
with an ammonia atmosphere, where the nickel powder was sintered to itself 
and to the nickel backup foil. The nickel powder used to make each sample 
is termed an irregular powder, irregular in that the nickel particles in the 
powder were more needle-like than spherical. The length-to -diameter ratio 
of these needle-like particles is typically 4 or less. 

In an irregular powder with its needle-like particles, the porosity can be 
varied over a fairly wide range from sample to sample. However, if a 
spherical powder is used; the porosity is restricted to fairly narrow limits. 
There are two regular patterns in which spheres of equal diameter can be 
packed together stably, cubic and rhombohedric. The porosities of these 
packings are 47. 84 and 25. 95 per cent respectively. By using irregular 
powders rather than spherical regular powders , the porosity of sintered 
powdered samples can be greater than the 47. 84 per cent that a spherical 
powder might impose. 

One of the sintered powder samples, M4, is shown in cross-section in Figure 
12 and in plan in Figure 13 . The solid structure at the bottom of the sample 
in Figure 12 is the nickel backup foil. The porous sintered nickel powder is 
0. 0873 inch thick and the nickel backup foil is 0. 0248 inch thick. 

The irregularity of the pore structure shown in Figures 12 and 13 is in con- 
trast to the relatively regular structure exhibited by the sintered screen 
samples. At some locations the powder seems to have melted in sintering 
and coalesced into large globules. 

The range of pore sizes is apparently greater in Sample M4 than for the 
sintered screens. As shown in Figure 13, the pores on the surface range 
from voids of dimensions comparable to the diameter of the smallest 
powder used, to large voids that could more aptly be termed craters. 

Figure 12 also shows that there is a significant number of large inter- 
connected pores of relatively low tortuosity running from the backup foil 
to the free surface of the porous material. Such a pore geometry is 
necessary for good boiling heat transfer characteristics, as will be shown 
in Section VIII. 
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Mag:21X 

Figure 12 Cross-Sectional View of Sintered Nickel Powder Sample M4 



Mag:21X 

Figure 13 Top View of Sintered Nickel Powder Sample M4 
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3. Sintered Fiber Materials 


The thirteen sintered metal fiber samples studied are listed in Table 3. 

TABLE 3 

Sintered Fiber Samples 


Mean Fiber 

Sample No. Material Diameter, Inches 


HI 

nickel 

0. 0004 

H2 

nickel 

0. 0004 

H3 

nickel 

0. 0004 

H4 

nickel 

0. 0004 

H5 

nickel 

0. 0006 

H6 

nickel 

0. 0006 

H7 

nickel 

0. 0006 

H8 

nickel 

0. 0006 

H9 

nickel 

0. 0006 

H10 

nickel 

0. 0006 

Hll 

430 stainless steel 

0. 0013 

H12 

430 stainless steel 

0. 0013 

H13 

430 stainless steel 

0. 0030 


As shown in the table, these sintered fiber samples were made from nickel 
or stainless steel. In each case the backup foil to which the porous material 
is sintered was the same material as the porous structure. 

The metallic fibers making up the porous material had a length-to-diameter 
ratio in the range of 10 to 1000. The fibers were felted to produce a randomly 
interlocked structure. The latter was pressed against a backup foil and then 
sintered in a eontrolled-atmosphere furnace to produce bonds at points where 
the fibers touch one another and the foil backing. 

A sintered nickel fiber sample, HI, is shown in cross-section in Figure 14 
and in plan in Figure 15. The nickel fiber structure is 0. 0825 inch thick and 
the nickel backup foil is 0. 0218 inch thick. 

Photomicrographs of a sample (H12) composed of stainless steel fibers of 
a much larger diameter are shown in Figures 16 and 17. In this case, the 
stainless steel fiber structure shown is 0. 0991 inch thick and the stainless 
steel backup foil is 0. 0218 inch thick. 
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Mag:21X 

Figure 14 Cross-Sectional View of Sintered Nickel Fiber Sample H 1 



Mag:21X 

Figure 15 Top View of Sintered Nickel Fiber Sample HI 
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Mag:21X 

Figure 16 Cross-Sectional View of Sintered Nickel Fiber Sample H12 



Mag:21X 

Figure 17 Top View of Sintered Nickel Fiber Sample H12 
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An examination of the photomicrographs in Figures 14 through 17 shows that 
there is an abundance of pores present in the fiber mat. In addition, there 
is a large degree of interconnectiveness associated with the pore structure. 
Because of their structure the sintered fiber metals offer the widest range 
of pore sizes and porosity of any of the porous materials that have been 
considered so far. 

B. Properties of Wicking Materials 

Knowledge of the porosity and the way in which the pore space is distributed 
through a porous solid is essential to any study of porous material. A des- 
cription of the quantities necessary to evaluate the porosity, the pore size dis- 
tribution, and the free-flow area ratio is presented in the following sections. 

1. Porosity 

The porosity, f, of a material is defined as the total pore volume of a sample 
divided by the total volume of the sample itself. Thus a porosity of zero 
per cent denotes a solid substance containing no pores or voids, while one 
of 100 per cent denotes a complete void. While the porosity accounts for 
all pores or voids in a sample, the effective porosity is based upon the 
interconnected (as opposed to separate or noninterconnected) pore space 
within the sample. 

The study of the flow of a fluid through a pore system in a material should 
be based on the effective porosity rather than the porosity. By examining 
Figures 10, 11, 14 and 15 it can be seen that the pores of both the sintered 
screens and the sintered fibers are for the most part interconnected. Thus 
the effective porosity of each of these materials can be taken as equal to the 
porosity. 

The powders used in making the sintered powder samples were composed 
of solid particles that contained no voids. However, in compacting the 
powder, isolated voids can form. This compacting was kept to a minimum 
during the fabrication of the sintered powder samples, represented by 
Figures 12 and 13, by compressing the powder before sintering as little as 
possible. Thus, in all samples considered in this study, the effective 
porosity will be taken as equal to the porosity, € . 

The porosity of each of the samples described in Section IV A was measured 
and is tabulated in Table 4. 
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TABLE 4 


Porosity of Porous Material Samples 


Sample No. Porosity, % 


Sample No. 



Ml 

Sintered 

M2 

Powder 

M3 


M4 

Samples ^ 

M5 


M6 


1 M7 

Sintered 

M8 

Screen 

M9 

Samples 

M10 


59.7 

65.8 

47.7 

54.0 

69.1 

69. 6 
62. 5 

67.9 

67.8 

67.6 


Sintered 

Fiber 

Samples 


HI 

H2 

H3 

H4 

H5 

H6 

H7 

H8 

H9 

H10 

Hll 

H12 

H13 


Porosity, % 


86.8 

82.5 
68.9 
62.8 

91.8 
88. 0 
82. 0 

82.8 

62.6 
70.9 
91.6 
80.8 
82.2 


In calculating the porosity values listed in Table 4 the basic definition of 
porosity was used in the form 


€ 


v _ W fp 
T T /m s 

V 


(55) 


where = total wick volume 

Wt = total weight of wick 

p s = density of the solid (8. 90 gm/cc for nickel and 7. 70 gm/cc 
for 430 stainless steel) 

The weight of each sample was determined by using a triple-beam balance, 
while the volume was determined by taking many measurements of the di- 
mensions of each sample with micrometers. In all measurements and com- 
putations, the fact that the samples were backed with a solid nickel or stain- 
less steel foil was taken into account. The maximum uncertainty of the per 
cent porosity (confidence interval based on odds of 20 to 1) for any of the 
sintered screen samples was calculated to be ±1. 6 while that of the sinter- 
ed powder and sintered fiber samples was ±1. 3. These measured porosity 
values not only provide a measure of the total pore space in a given sample, 
but also tell something about the structure of the porous material. For in- 
stance, if the sintered powder samples were composed of perfectly spher- 
ical powders, the maximum porosity that each sample could have would be 
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47. 84 per cent. However, the substantially higher values of porosity listed 
in Table 4 show that the combined effect of irregular powders and loose pack- 
ing leads to a structure that is less dense than an equivalent spherical packing. 

The maximum porosity that a perfectly stacked screen structure can assume 
was calculated to be about 61 per cent. This configuration was termed a 
square stacking. The measured porosity values of Table 4, ranging from 
62. 5 to 67. 9 per cent, indicate that the screen samples in this study probably 
approximate the square stacking. Although Figure 10 shows that many of 
the layers follow the square configuration, there are many layers that are 
much more closely stacked. The higher porosity values that were measured 
indicate that some of the screen layers may be imperfectly sintered, thus 
producing a greater void volume in the sample. 

The sijliered fibers have a much higher porosity than the other two classes 
of porous material, due to their structure. In one case the randomly inter- 
locked fiber structure produced a porosity of 91. 8 per cent. This is about 
the practical upper limit on rigid sintered fiber materials, for this sample 
(H5) tended to deform easily. The porosity of any of these samples, be 
they sintered screen, powder, or fiber, could be decreased by rolling and 
compressing the sample. 

2. Pore Size Distribution 


The manner in which the pore space is distributed in a porous material is 
given by the pore size distribution function. This function, a(r), is defined 
by the equation 


dV = - a(r)dr (56) 

where dV is the volume of the pores with radius between r and r + dr. For 
purposes of the following discussion, V itself will be defined as the volume 
of all pores with radius equal to or greater than r. 

The porous materials described in Section IV A have pore diameters falling 
in the range of 0. 1 to 1000 microns. For this range of pore diameters the 
best method available for determining the pore size distribution function is 
to use a mercury intrusion porosimeter. This instrument employs the 
injection of non -wetting mercury into the pore system of a sample, and 
recording the pressure necessary to accomplish the injection. By measuring 
the volume of the mercury in the porous material as a function of pressure, 
the pore size distribution can be calculated. An explanation of the opera- 
tion of the device is presented in the following paragraphs. 


When a liquid is introduced into the pore system of a porous material, the 
liquid will form menisci whose curvature is determined by the shape and 
size of the pores in the material. If the liquid is a non-wetting fluid such 
as mercury, a positive pressure will be necessary to force the liquid into* 
the pore system. For a given size of pore, this positive pressure is given 
by the Laplace - Young equation 



(57) 


where AP is the pressure drop across the liquid-vapor interface, <r is the 
liquid-vapor surface tension, and Rq and R 2 are the two principle radii of 
curvature necessary to describe the liquid-vapor interface. If the pore 
system has been evacuated prior to the introduction of the liquid and if the 
vapor pressure of the liquid is small (e. g. , mercury at room temperliire), 
then 


P ^ AP 


(58) 


where P is the pressure of the liquid. 


It is assumed here that Ri and R 2 completely describe the entire liquid- 
vapor interface in the pore. Up to now the terms pore radius and pore 
diameter have been used without precise definition. An examination of the 
photomicrographs, Figures 10 to 17» shows that it would be quite unrealistic 
to assume that the pores are cylindrical in the materials under study. Their 
pore shape and structure are quite irregular. However an equivalent pore 
radius can be defined as 

- 2 ( cross ~ sec ^ Qna ^ urea of pore 
e \ perimeter of pore 

In fluid mechanics terminology this equivalent pore radius can also be called 
a hydraulic radius*. As r e is defined, Equation (59) reduces to the radius of 
a cylinder for the case of a cylindrical pore. Hereafter the term pore radius 
and r will be taken as equivalent to the definition given by Equation (59). 



Using Equation (58), a force balance on the meniscus in the pore yields 

(F a ) (perimeter of pore) = (P) (cross-sectional area of pore) (60) 

where F a is the surface tension force per unit length of interfacial boundary. 
The latter is given by 


F a = a cos <t> (61) 

where a is the surface tension and <f> is the contact angle. 
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Combining Equations (59), (60) and (61) yields 

Pr = 2<rcos0 (62) 

For constant values of surface tension and contact angle, this equation 
directly relates r, the pore radius, to P, the measured pressure necessary 
to force the non-wetting liquid (in this case mercury) into the pore. It 
should be noted that the expression is independent of the measured volume 
of liquid, V. If 0 and c are constant, Equation (62) can be differentiated 
to obtain 



Thus to evaluate a( r) for a given porous material, both P and V must be 
measured. This is done with the mercury intrusion porosimeter by meas- 
uring the pressure, P, necessary to force a given volume ,V, of mercury into 
the evacuated sample. Knowing this a(r) can be calculated from Equation (64) 
and r can be found from Equation (62), where a is taken to be 473 dynes/cm 
for mercury at 70°F and the contact angle is 130 degrees. The latter value has 
been found to be valid for most mercury-metal systems. The derivative dV/dP 
is calculated by taking the slope of a plot of P vs V. 


PAGE NO. 41 



Figure 18 Pore Size Distribution Curve 


A typical pore size distribution curve is shown in Figure 18, where r is 
defined as the mean pore radius, i. e. , the pore size that occupies more of 
the total pore volume than any other pore size. It should also be noted that 
as defined, a(r) is dependent on the actual size of the sample. 

The V-P data to calculate a(r) was obtained for twelve of the 
most promising porous material samples described in Section IV A. These 
were samples HI, H2, H3, H6, H10, Hll, H12, H13, M2, M4, M6, and M8. 
At least two pieces of each porous material sample were tested in a mercury 
intrusion porosimeter made by the American Instrument Company. Each 
sample was placed in the porosimeter container which was evacuated before 
the mercury was introduced into the sample. All samples were cleaned 
using the procedure outlined in Section V C. The most complete V-P curve 
for each of twelve samples is shown in Figures 19 through 30. These plots 
show volume as a function of pressure and volume as a function of pore 
diameter. From Equation (62) the pore diameter is given by 

= 4 g- cos <p_ , 6S . 
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where 


D = diameter of pore, i. e. , twice the pore radius given by 

Equation (59) 

P = pressure required to force mercury into a pore of 
diameter D 

<r = surface tension of mercury = 473 dynes/cm at 70 °F 

<t> = contact angle between mercury and solid material = 

130 degrees 


Figures 19 through 30 present essentially the same information that a plot 
of the pore size distribution would show, for they illustrate the measurable 
range of pore sizes in each sample and the volume occupied by each pore 
size. To replace the mean pore radius r in Figure 18, the median pore 
diameter has been chosen to characterize the pore space of each sample. 

The median pore diameter D me d is defined as the pore size at which 50 per 
cent of the connected pore space is filled with mercury. 

Fairly complete pore spectra were obtained for samples HI, H2, H3, H6 and 
M10, with pore diameters falling between 0. 5 and 100 microns (Figures 19 
to 23). These samples were four of the sintered fiber samples of finest 
fiber and one sintered powder sample. The range of pore sizes that occu- 
pies the greatest portion of the total pore volume occurs in the region of the 
maximum slope of the V-P curve. For sample H3 (Figure 21) this portion 
of the total volume occurs for pore sizes between 20 and 30 microns. For 
this sample the median pore diameter Dmed can be calculated by noting that 
at a pore volume of 0. 08 cc (equal to 50 per cent of the connected pore 
volume) the pressure as read from the V-P curve of Figure 21 is about 6.6 
psia. Using Equation (65), or using the pore diameter scale in Figure 21, 
D m ed * s f° un d to be 26. 5 microns. Figure 21 also shows that the pore 
diameters ranging from 10 to 34 microns contain 90 per cent of the connected 
pore volume of the sample. Similar conclusions can be drawn from Figures 
19, 20, 22 and 23. 

Because of the construction of the mercury intrusion porosimeter, pressures 
less than 2. 5 psia could not be read. This pressure, using Equation (65), 
corresponds to a pore size of 70 microns. Thus, if a sample had pore 
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structure in which most of the pore space was contained in pores of diameter 
greater than 70 microns, a complete pore spectrum could not be measured 
by this method. 

Such was the case in the remainder of twelve selected materials, samples 
Hll, H12, H13, M2, M4, M6 and M8. They constitute the more porous 
sintered fiber samples, two sintered powder samples and one sintered 
screen sample. The V-P curves for these samples are shown in Figures 
24 to 30. In all of these cases the median pore diameter is greater than 70 
microns. However, because of the shape of the V-P curves shown in 
Figures 24 to 30, the value of the lower bound estimate on median pore size 
can be increased. For example, sample M8 can be interpreted to have a 
median pore size greater than 100 microns. 

TABLE 5 

Pore Space Distribution of Twelve Selected Porous Materials 
Measured Median 

Pore Diameter, Pore Diameter Range of 

Sample No. microns 90% of Pore Space, microns 


M2 

>70 

- 

M4 

- 80 

- 

M6 

> 80 

- 

M8 

>100 


M10 

~ 50 

43.4-57.8 

HI 

43.8 

19.4-54.2 

H2 

35. 0 

13.3-48.9 

H3 

26. 5 

10.3-34.0 

H6 

40.7 

15. 5-59.3 

Hll 

>80 

- 

H12 

>70 

- 

H13 

>80 

— 


The pore space distribution data presented by Figures 19 to 30 is summarized 
in Table 5. The median pore diameter is listed for each sample. Also 
tabulated for five of the twelve samples is the pore diameter range that 
represents 90 per cent of the pore volume in any given sample. In all cases 
the term pore diameter is defined as twice the radius defined by Equation (59). 

3. Free-Flow Area Ratio 

The free-flow area ratio is defined as the area available to flow the fluid 
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Figure 19 Mercury Intrusion Porosimeter Curve for Sample HI 
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Figure 20 Mercury Intrusion Porosimeter Curve for Sample H2 
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Figure 21 Mercury Intrusion Porosimeter Curve for Sample H3 
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Figure 22 Mercury Intrusion Porosimeter Curve for Sample H6 


PAGE NO. 48 



PENETRATION VOLUME, 



PORE DIAMETER, D~ MICRONS 


IOO 50 20 10 5 2 1.0 0.5 0.2 



ABSOLUTE PRESSURE, P~ PSI 


Figure 23 Mercury Intrusion Porosimeter Curve for Sample M10 
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Figure 24 Mercury Intrusion Porosimeter Curve for Sample Hll 
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Figure 25 Mercury Intrusion Porosimeter Curve for Sample H12 
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Figure 26 Mercury Intrusion Porosimeter Curve for Sample H13 
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Figure 27 Mercury Intrusion Porosimeter Curve for Sample M2 



PAGE NO. 53 



PENETRATION VOLUME, V~cc 



ABSOLUTE PRESSURE, P~PSl 


Figure 28 Mercury Intrusion Porosimeter Curve for Sample M4 
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Figure 29 Mercury Intrusion Porosimeter Curve for Sample M6 
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Mercury Intrusion Porosimeter Curve for Sample M8 
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at any given cross-section of the porous material, divided by the total 
cross-sectional area. In accordance with this definition, the free-flow area 
ratio is a measure of the actual cross-sectional area available to the flow of 
a fluid in any given porous material. 


i.n expression for F^r that is based solely on considerations of isotropy can 
be derived as follows. In a given porous material, a body of any shape can 
be considered so that its smallest dimension is much greater than the largest 
pore size. The body will be considered to be xx units long in some chosen 
x-direction and to have a cross-sectional area normal to the x-direction of 
Ax, where Ax is some known function of x. 


By definition, the free-flow area ratio at a given x for the body of porous 
material is given by 


Alt 



( 66 ) 


or 


= F A 
Alt T 


(67) 


where A^ is the total cross-sectional area at x, F^r is the free-flow area 
ratio, and Ap is the total pore area contained in A T . 


The porosity of the body is given by 



( 68 ) 


where e is the porosity, Vx is the total volume of the body of porous material, 
and Vp is the total volume of the pore space in the body. The total volume can 
be calculated by integrating in the x-direction to get 



Ax dx 


(69) 


Similarly, by using Equation (67) for Ap, the pore volume is given by the 
integral 



( 70 ) 
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Combining Equations (68), (69), and (70), it is seen that 


r x l 

J o Far A t dx 



( 71 ) 


As stated above, A T is some known function of x. Now, if it is assumed that 
the porous material is truly isotropic and that the pore space is randomly dis- 
tributed, the value of the free- flow area ratio as given by Equation (66) at any 
given x will be a constant, if the dimensions of A^ are large compared to the 
largest pore sizes. Thus, since F^j^ is not a function of x. Equation (69) 
reduces to 


f AR - f ( 72 ) 

subject to the conditions of isotropy and a randomly distributed pore space. 
Hence, for an isotropic porous material with a randomly distributed pore space, 
the free-flow area ratio is equal to the porosity. This result agrees with a 
similar conclusion given by CarmanlO-. 

The photomicrographs of Figures 12 to 17 indicate that the assumptions of iso- 
tropy and a randomly distributed pore space apply reasonably well to sintered 
powder materials and sintered fiber materials. However, the structural 
regularity exhibited by the sintered screen sample shown in Figures 10 and 11 
indicates that Equation (72) is at best an approximation to a mean value of 
FAR f° r these porous materials. The interlocking of the screen layers dis- 
cussed in Section IV A will help both assumptions. 

The values of the free-flow area ratio are useful in computing actual fluid flow 
velocities, Reynolds numbers, and an equivalent flow pore diameter. 

C. Porous Material Machining 

One of the problems encountered when dealing with porous materials was that 
of cutting or machining a sample without disturbing the porous structure in the 
region adjacent to the cutting plane. A minimum of disruption is necessary to 
insure that there are no end effects at the edges of the samples used for the 
wick permeability and boiling studies. 

Normal cutting techniques such as band- sawing, shearing, grinding and cutting 
with a torch were found to both deform the porous material sample itself and 
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smear the cut surface. Such a surface is shown in Figure 31, a photomicro- 
graph of a cross-section sample H4, a sintered nickel fiber material. The 
machined surface shown was produced by cutting the sample with a rotating 
rubber wheel charged with diamond dust. The nickel backup foil (0. 0218 inch 
thick) is at the bottom of the cross-section, with porous structure (0. 0982 inch 
thick) above it. It can be seen that almost the entire pore space has been 
smeared over and sealed off by the machining process. Such a surface would 
exhibit a relatively high resistance to flow during a permeability test. 

The porous structure can be preserved during cutting by filling the sample with 
some supporting material such as a wax or a salt. The filler material can be 
removed after the machining operation by melting or leaching. However, it is 
difficult to determine whether the filler is completely removed or that it has 
not changed the wetting characteristics of the porous material. 

Figure 32 shows a photomicrograph of the same sample H4 where the material 
was cut by electro-discharge machining (EDM). EDM is a process of metal 
removal brought about by an electrical spark jumping a gap between the work- 
piece and an electrode. A dielectric liquid flows between the electrode and the 
workpiece, both to cool the work and to remove the spark-eroded particles as 
they come off the workpiece. 

As seen by the cross-sectional view of Sample H4 in Figure 32 , this machining 
method tends to produce a rough edge. The electrical discharges tend to erode 
the solid backup foil unevenly, producing a pitted surface. When examined under 
a microscope the porous material near the cut surface can be seen to be dis- 
rupted and opened up just the reverse of the situation shown in Figure 31. 

The cutting process finally chosen was electrochemical machining (ECM). A 
cross-sectional view of Sample H4 cut in this manner is shown in Figure 33. 
Electrochemical machining is essentially metal removal by controlled electro- 
plating in reverse. The workpiece and a shaped tool form a pair of electrodes. 

An electrolyte (usually an ordinary salt solution) is pumped between the work- 
piece and the tool, and a voltage is applied to them. The material of the work- 
piece is then removed by electrolysis, duplicating the shape of the tool. The 
electrolyte sweeps away the metal ions so that they cannot deposit out on the 
tool, which acts like an anode. This cutting process occurs at the temperature 
of the electrolyte so that there are no excessive temperatures at the cut surface. 
Also, because of the electrolysis process, all chemical reactions take place 
right at or near the cut surface rather than in the matrix of the porous material. 
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Mag:21X 

Figure 31 Cross-Sectional View of Sintered Nickel Fiber Sample H4 Cut 
by Grinding Method 



Mag:21X 

Figure 32 Cross-Sectional View of Sintered Nickel Fiber Sample H4 Cut 
by Electrical Discharge Method 
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Figure 33 shows that both the porous material and the backup foil of Sample H4 
were cut cleanly and with a minimum of disruption. As will be shown in Section 
VI D , samples cut by the ECM method show no observable end effects. The 
machined surfaces shown in Figures 10, 12 , 14 and 16 were all produced by 
the electrochemical machining process. 



Mag:21X 

Figure 33 Cross-Sectional View of Sintered Nickel Fiber Sample H4 Cut 
by Electrochemical Machining Method 

D. Test Fluids 


As shown by Equation (45), the liquid properties enter into the operation of the 
vapor- chamber fin as the dimensional parameter N, defined as 

N = PL hVL<r (73) 

where Pl is the liquid density, hyL is the latent heat of vaporization, a is the 
liquid-vapor surface tension, and Pl is the liquid viscosity. 

To maximize the heat transport properties of the fin, N should be as large as 
possible. Also, the test liquid should have a vapor pressure significantly higher 
than the vapor pressure drop in the vapor-chamber fin. In this way the tempera- 
ture of the vapor-chamber fin can be maintained nearly uniform, resulting in 
an efficient device. Also, the chamber pressure should not be too high, or the 
fin will become a pressure vessel with prohibitively thick walls. 
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Values of N for a number of liquids are listed in Table 6. 


TABLE 6 



Properties of Various Liquids 
Normal 

Boiling Vapor Pressure, psia N*, 

Btu hr lbf/ft 3 

lbm 

Liquid 

Point, °F 

at 70°F 

at 1300°F 

at 70°F at 1300°F 

at B. P. 

Metal 

bismuth 

2875 

_ 

<10" 2 

_ 

2097 

2360 

lead 

3169 

- 

<10” 2 

- 

2007 

- 

sodium 

1630 

- 

2. 0 

- 

1983 

2016 

mercury 

674 

- 

760 

- 

- 

1253 

potassium 

1395 

— 

9. 0 

— 

406 

411 

Non-metal 

water 

212 

0. 36 

_ 

138 

_ 

336 

methyl alcohol 

150 

- 

- 

25. 1 

- 

- 

acetone 

136 

3.77 

- 

17.9 

- 

- 

ethyl alcohol 

173 

- 

- 

9.9 

- 

19.2 

Dowtherm E 

352 

<5.0 

- 

9. 8 

- 

- 

carbon 

tetrachloride 

- 

1. 8 

- 

7.3 

- 

- 

Freon 112 

199 

- 

- 

- 

- 

5.7 

Freon 12 

-21. 6 

84 

- 

5.2 

- 

- 

Freon 113 

118 

5. 5 

- 

4. 8 

- 

5.2 

Freon 114B2 

117 

5.5 

- 

4.2 

- 

4.4 

Dowtherm A 

496 

<0.2 

- 

2. 9 

- 

- 

FC-75 

216 

0.4 

- 

1. 15 

- 

1.21 

sulfur 

832 

- 

245 

- 

- 

0. 17 


* approximate values 

The values listed in this table show water to have one of the highest values of 
N for conventional liquids, in addition to having a favorable vapor pressure - 
temperature relationship for testing. Hence water was selected as the primary 
test fluid. 

To evaluate the effect of the variation of N on the heat transfer properties of 
the fin, one other liquid was selected. The second fluid must satisfy the condi- 
tions given above and have a value of N significantly lower than that of water, 
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although not low enough to drastically reduce the heat input to the fin, see 
Equation (49), and thereby cause the experimental uncertainties to become too 
high. Also the fluid should have a fairly complete tabulation of properties and 
be compatible with the wicking material. Because it was to be tested in the at- 
mosphere during the wicking and boiling tests it must be relatively nonreactive 
and nontoxic. 

Of the liquids listed in Table 6, Freon 113 seems to best satisfy these require- 
ments and thus was used as the second fluid to evaluate the effect of N on the 
operation of the vapor- chamber fin. 

E. Summary 

A total of 23 porous materials was selected, representing three classes of 
porous materials: sintered metal screens, sintered metal powders, and sint- 
ered metal fibers. Each sample was composed of either nickel or 430 stainless 
steel. 

The porosity of all the samples was measured and was found to vary between 
47. 7 and 91. 8 per cent. 

A mercury intrusion porosimeter was used to measure the pore size distribution 
function of twelve of the twenty-three porous material samples. A complete 
pore spectrum was found for five of the twelve samples. The median pore 
diameter for these five ranged from 26. 5 to 50. 0 microns. 

The free-flow area ratio, useful for computing Reynolds number and effective 
flow velocities in a given porous material, was found to be equal to the porosity. 
Thus the free-flow area ratio ranged from 47. 7 to 91. 8 per cent for the twenty- 
three samples. 

A study of the machining of porous materials was made. It was found that the 
best method available was that of electrochemical machining and therefore this 
method was used in fabricating all samples used in subsequent testing. 

Water was selected as the primary test fluid since analysis indicated that it is 
one of the best conventional fluids for use in a vapor-chamber fin. Freon 113 
was selected as the secondary test fluid since its properties fall in the desired 
range and are significantly different from those of water. 
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A. Objective 


V. WICKING RISE TESTS 


The purpose of the wicking rise tests was to determine the maximum height to 
which a liquid will rise in a given vertical wicking material sample and there- 
by obtain a measure of the minimum liquid-vapor interfacial radius of curvature 
that the material will support. As shown in the derivation of the fin model in 
Section III D , the value of this curvature, Rmin. is needed to evaluate the maxi- 
mum pumping force produced by the wicking material. This maximum pumping 
force, expressed as a differential pressure is 

2 <r 

AP = — (74) 

“mm 

Thus by evaluating Rmin for a given wicking material, the maximum pumping 
force of the wick is evaluated. The frictional forces associated with the liquid 
flow caused by the pumping force are evaluated in the permeability tests des- 
cribed in Section VI. 

The wicking rise test is a simple method of determining the effective minimum 
radius of curvature that a wicking material can support. As an explanation, 
the wick and liquid reservoir shown in Figure 34 can be considered. 



Figure 34 Equilibrium Condition for Three Selected Pores in Saturated Wick 
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In the porous material shown in Figure 34, only three connected pores are pic- 
tured out of the many that actually contain liquid. It is assumed that the material 
of this wick is one of those listed in Section IV A. Thus, due to the properties 
of isotropy and homogeneity, at any given height z in the wick, all pore sizes 
from the smallest to the largest will be represented. 

In Figure 34 the liquid has risen to an equilibrium height in each of the pores 
and it is observed that the liquid front in the wick has risen in some pores to 
a maximum height Jtm represented by pore C. It has been noticed that in a 
wicking rise test, the degree of wick saturation decreases with increasing z. 

This is probably due to the smaller fraction of pores (i. e. those with sufficiently 
small radii) being able to overcome the increased liquid head as z increases. 

In other words, as z increases, more ’and more of the larger pores are unable 
to support the liquid head. 

By employing the laws of hydrostatics in the equilibrium situation pictured, the 
pressure in the liquid of pore C at point 1 just under the meniscus is 

Pl_,S n 

Pl = Pa - — < 75 > 


where P a is the pressure on the liquid reservoir at z = 0. Equation (75) holds 
no matter what the shape and size of pore C is, or how many connections it has 
with other pores and the free surface of the wick. 

If the very small vapor head between z = 0 and point 2 just above the meniscus 
is neglected, then 


p 2 = 


The pressure difference across the meniscus in pore C is then 


(76) 


P 2 " P l = 




S r 


m 


(77) 


But the Laplace- Young equation (see Section III D ) gives the pressure differ- 
ence across the meniscus as 


P 


2 



2 a_ 
R c 


(78) 


where it is assumed that pore C can be characterized by one value of the radius 
of interfacial curvature, R c . 
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By equating Equations (77) and (78) it is seen that 


p = Q * — (7Q\ 

R C p gft ( 79 > 

JLj m 

Since is a maximum for the sample and all pore sizes are present at any 
given z 


R c R min 


(80) 


where Rmin is the minimum radius of curvature that the wicking material can 
support. Thus the wicking tests are a means of determining R min . 


The above assumes that the liquid wets the porous material of the 
wick, i. e. , has a contact angle 0 < <t> < 90°. The wicking test is not a mea- 
sure of the minimum pore size of the wicking material, because a definite 
value of the contact angle would have to be known to calculate this. It merely 
permits an evaluation of the maximum capillary forces that a given liquid- 
wicking material combination will produce. 


Thus two wicking materials with exactly the same pore structure but made of 
different materials would probably have different equilibrium heights with a 
given liquid, due to differences in their wetting characteristics. Also if a 
given wick is considered with two different liquids, the minimum radius of 
curvature, Rmin> °f each wick- liquid combination could be different due to 
differences in the wetting characteristics. 


The result of the wicking rise tests is an equilibrium height yielding a measure 
of the maximum capillary pumping force for the wick- liquid combination tested. 
The pumping forces are functions of both the wetting characteristics of the wick- 
liquid combination and the pore size and structure of the wick. 


B. Description of Wicking Rise Apparatus 

The test apparatus consists of glass tubes with a 60 mm inside diameter which 
are attached in a vertical position to a metal framework as shown in Figure 35. 
Each tube encloses a wicking sample and is open at the bottom with a rubber 
stopper on the top. The stopper has two openings in it to receive a valve assem- 
bly, used in venting the tube before the start of a test, and a fitting which allows 
the sample height to be adjusted. The sample is suspended in the glass tube by 
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Figure 35 Wicking Apparatus with Samples Installed 
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a stainless steel rod and clip that passes through the fitting. Also attached to 
each rod is a thermometer used in recording the temperature in the tubes. 

The glass tubes are suspended in a Plexiglas reservoir which has a Plexiglas 
top fitted with soft Neoprene seals that fit around the outside of each glass tube. 
The seals and the Plexiglas top prevent evaporation and contamination of the 
liquid contained in the reservoir. The Plexiglas reservoirs are large enough 
so that the level will remain nearly constant while tests are in progress. The 
reservoirs are filled with either water or Freon 113, the two liquids under 
investigation. 

The bottom of each tube is submerged below the surface of the liquid in the re- 
servoir, thereby sealing the tube and sample from the atmosphere. The air 
in the tube is saturated, preventing any net evaporation from the wick during 
and after the time the liquid has risen to the equilibrium height. 


The results of the wicking tests have shown that it can take a long time, i. e. , 
days, for the liquid front in a wicking sample to reach an equilibrium height. 
Some of the denser porous metal samples did not reach an equilibrium height 
in 20 days. Similar results have been obtained by other investigators^. Gen- 
erally, after the first few hours of a wicking test, the rate of rise of the liquid 
front becomes very small and increasingly difficult to see. Because of the 
small amount of liquid present at the higher levels in the sample, it is quite 
important to enclose the wick in a saturated- vapor environment to prevent mass 
transfer. Failure to eliminate the latter would cause a lower than actual value 
of j? m , the equilibrium height, to be recorded. During tests each apparatus 
was completely enclosed in a multiwalled cardboard enclosure to keep the entire 
apparatus at a near constant temperature and to minimize temperature gradients 
along the length of the wicking samples. As shown in Figure 36, each enclosure 
had a Plexiglas front for viewing the samples during a test. 

During earlier tests it was found that the temperature gradient along each tube 
could amount to as much as 4°F. This not only caused errors in estimating 
surface tension but also caused condensation to form on the inside of each tube, 
making the reading of the liquid level in each sample difficult. To remedy this, 
thermocouples were attached to the outside of three of the tubes, at the top, 
middle, and bottom. Thermocouples were also placed in each of the two liquid 
reservoirs. After the enclosure described above was installed, the tempera- 
ture gradient along each tube decreased to less than 1°F, and condensate no 
longer formed in the tubes. 
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Figure 36 Wicking Apparatus with Insulating Enclosure 
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Due to the surface finish and the opacity of the materials tested, it was found 
that it was not possible to visually observe the liquid front in a given wick. 

Early in the experimental investigations various methods of measuring the posi- 
tion of the liquid front were tried, the best one being that of attaching litmus 
paper to the wick and noting when the paper was first wet. Other methods tried 
were the use of dyes and fluorescent liquids irradiated with ultraviolet light. 

The litmus paper method was eventually refined by adding a small amount of 
potassium hydroxide to the water to increase the response of the litmus paper. 
Tests were made to find the smallest amount of basic material that could be 
added to the water and still cause pink litmus paper attached to a wick to turn 
to a vivid blue. It was found that a solution 0. 03 per cent by weight of KOH 
was adequate for this result. This wicking liquid can be treated as being essen- 
tially pure water, due to the very small amount of KOH used. Figures 37 and 
38 show that the interfacial properties of aqueous solution of 0. 03 per cent of 
KOH differ insignificantly from those of pure water. 

The litmus paper in the form of thin strips was held against the surface of the 
porous sample by a Plexiglas ladder. The ladder and the sample were held 
together by rubber bands placed around the assembly at various points. A 
picture of a wick sample both with and without the ladder is shown in Figure 39. 
The Plexiglas ladder shown in the sketch below, was made by machining rectang- 
ular grooves into a 1/ 8- inch thick sheet of acrylic plastic, resulting in lands 
about 0. 08- inch wide for holding the strips of litmus paper. The first ladders 
had a land pitch of 0. 5 inch, but later the pitch was set at 0. 3 inch. On a 
sample 20 inches long there were 66 strips of litmus paper approximately 0. 1 
inch wide, one under each land. The grooves in each Plexiglas ladder were 
0. 08 inch deep, making the ladder very flexible. When rubber bands were 
wrapped around it and the sample, the ladder conformed to the local irregular- 
ities of the porous sample, thereby pressing the litmus strips against the sam- 
ple surface for good contact. 
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KOH CONCENTRATION 

PER CENT BY WEIGHT OF KOH 

Figure 38 Equilibrium Contact Angle vs Aqueous KOH Concentration for Two 
Nickel Surfaces 
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There are two drawbacks to the use of the litmus paper. One is that the litmus 
paper does not respond instantaneously to the position of the actual liquid front 
when the front first reaches the level of a given strip. It takes a finite amount 
of time to wet the paper and to react with the pH indicators in the litmus paper. 

The measured rate-of-rise curve is then lower than the rate-of-rise curve of the 
actual liquid front in the sample. However, for elapsed times of an hour or 
longer, this difference was found to be negligible. 

The other drawback is the fact that each strip is separated by a finite distance 
of 0. 22 inch. Thus the finest division for measuring the height of the liquid 
front is this spacing. It was felt that any closer spacing would introduce too 
much litmus paper, which in itself is an excellent wick and might influence test 
results. However, with moderate wicking heights of 6 inches, this error of 
-0. 22 inch +0. 08 inch amounts to an uncertainty of only -3. 7 per cent +1. 3 per cent 
in measured values of maximum wicking rise, ^ m . 

Two wicking rise stands were built in order to test up to eleven wicking samples 
at a time. The samples varied in length from about 4 to 16 inches depending on 
the individual wicking capability. Each wick sample was sintered to a backup 
plate to simulate the actual construction in a vapor-chamber fin application. 

C. Wicking Rise Test Procedure 

The procedure involved in observing the liquid front in a wicking sample consisted 
of placing strips of pink litmus paper every 0. 30 inch along the vertical length of 
the sample. In the tests involving water the wicking liquid as explained in the last 
section was actually an aqueous solution of 0. 03 per cent potassium hydroxide, 
essentially having the same surface chemistry characteristics as pure water. 

The liquid front which is difficult to visually observe in the sample itself, was 
located by noting when each successive litmus paper started to turn blue. 

The same technique was utilized during the wicking tests in which Freon 113 
was used, except that blue litmus paper was substituted for the pink. Nothing 
was added to the Freon, for it was found that the blue litmus paper turned a 
very visible darker blue immediately upon coming in contact with Freon 113. 

Before any wicking test was begun, each sample was subjected to the following 
cleaning procedure: 

1) Washing the sample in a vapor degreaser and immediately rinsing 
the sample in distilled water before the condensed vapor 
re- evaporated. 


PAGE NO. 73 



2) Bathing in a reagent-grade acetone bath, 

3) Submerging the sample immediately in distilled water, followed 
by several distilled water rinsings, and 

4) Air-drying the sample in a clean oven at 500°F for at least 2 hours. 

This procedure was followed for each sample to try to insure that the wetting 
characteristics of each sample would stay constant. After a sample was re- 
moved and cooled in the atmosphere, it was stored in a clean plastic bag until 
tested. 

At the start of a wicking test, the sample with attached litmus paper was placed 
into one of the tubes shown in Figure 35, with the end of the sample above the 
reservoir level. Time was permitted for the air in the tube to come into equi- 
librium with the liquid so that there was no net mass transfer during the actual 
wicking test. The lower end of the vertical wicking sample was lowered into the 
liquid by means of the fitting at the top of the tube , and a stopwatch started. The 
time, temperature, and level were recorded every time a litmus paper changed 
color, indicating passage of the liquid front. When the rate of rise slowed to 
less than 0. 3 inch (one litmus paper) per hour, readings were taken every hour. 
When the rate of rise of the liquid front was less than 0. 3 inch for every two 
hours, readings were taken at two-hour intervals to the end of a test. The 
wicking run on a sample was terminated when either, 1) the liquid front reached 
the sample top, 2) a period of two weeks elapsed, or 3) the liquid front failed 
to rise 0. 3 inch (one litmus paper) in 3 days. The height at which the latter 
occurs is defined as the equilibrium height, X m . 

D. Test Results 

At least one wicking rise test was run on each of the samples described in 
Section IV A. In some cases as many as three wicking tests in water and one 
in Freon 113 were run on a sample. 

The results of the wicking tests in water are shown in Figures 40 to 61. 
which present the variation of liquid height with time. The only sample listed 
in Section IV A that is not represented in these figures is Sample H5. This 
sample was excluded after it was found that it had not been completely oven- 
dried after the cleaning process, which resulted in erroneous data for the only 
wicking test run on H5. In all cases the curves shown in these figures repre- 
sent the particular wicking test that resulted in the largest value of the equilibrium 
height. 
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Figure 40 Water Wicking Curve for Sintered Nickel Power Sample Ml 
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Figure 41 Water Wicklng Curve for Sintered Nickel Power Sample M2 
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Figure 42 Water Wicking Curve for Sintered Nickel Power Sample M3 
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Figure 43 Water Wickiug Curve for Sintered Nickel Power Sample M4 
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Figure 44 Water Wieking Curve for Sintered Nickel Power Sample M5 
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Figure 45 Water Wieking Curve for Sintered Nickel Power Sample M6 
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Figure 46 Water Wieking Curve for Sintered Nickel Screen Sample M7 
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Figure 47 Water Wieking Curve for Sintered Nickel Screen Sample M8 
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Figure 48 Water Wicking Curve for Sintered Nickel Screen Sample M9 



Figure 49 Water Wicking Curve for Sintered Nickel Screen Sample M10 
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Figure 50 Water Wicking Curve for Sintered Nickel Fiber Sample HI 
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Figure 51 Water Wicking Curve for Sintered Nickel Fiber Sample H2 
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re 52 Water Wicking Curve for Sintered Nickel Fiber Sample H3 
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Figure 53 Water Wicking Curve for Sintered Nickel Fiber Sample H4 
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Figure 54 Water Wicking Curve for Sintered Nickel Fiber Sample H6 
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Figure 55 Water Wicking Curve for Sintered Nickel Fiber Sample H7 
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Figure 56 Water Wicking Curve for Sintered Nickel Fiber Sample H8 
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Figure 57 Water Wicking Curve for Sintered Nickel Sample H9 
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Figure 58 Water Wicking Curve for Sintered Nickel Fiber Sample H10 



Figure 59 Water Wicking Curve for Sintered Stainless Steel Fiber Sample Hll 
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Representative wicking samples were chosen from those that had the largest 
value of .£ m with water. These samples, Ml, M4, M6, HI, H2, H3 and H 6, 
were then recleaned and prepared for wicking tests with Freon 113. The rate 
of rise results for four of these samples are shown in Figures 62 to 65. Although 
the equilibrium height was obtained for the other three samples, their rate of 
rise curves are not shown since the liquid reservoir height varied during the 
test. This variation was due to evaporation of the Freon 113 over weekends. 

Because of the relatively long time required for the liquid front to reach equilib- 
rium in a sample, all of the plots shown in Figures 40 to 65 have logarithmic 
time scales. As shown in these figures this equilibrium time can exceed 10, 000 
minutes. 

In a number of the tests the samples proved to be too short, as the liquid front 
rose to the top. This occurred with samples M3, M4, M8, HI, H2, H3, H4, 

H6, H7, H8, H9 and H10 in the tests involving water. Unfortunately, wicks 
longer than 16 inches could not be obtained from all of the wick manufacturers. 

In tests where the liquid front did not progress to the top of the samples (samples 
Ml, M2, M5, M6, M7, M9, M10, Hll, H12, H13 with water and all the tests 
involving Freon 113), the height of the first litmus paper above the liquid front 
and the time at which the test was terminated are indicated by a dotted circle. 

The solid circle just below the dotted circle indicates the last litmus paper to 
be wetted. This solid circle represents the equilibrium height of the sample 
except for the test on Sample M4 in Freon 113 (see Figure 63). This test was 
discontinued before the equilibrium height (as defined in Section V C) was 
reached. In the test on Sample Hll in water (see Figure 59) the observed height 
of liquid fluctuated with time near the final equilibrium height. This was per- 
haps due to a combination of difficulty in observing the height and to small tem- 
perature variations causing non-equilibrium conditions to exist. 

No attempt was made to fit an analytical model to any of the rate of rise curves 
shown in Figures 40 to 65. This is because these curves do not present the true 
rate of rise of the liquid front in the sample, but are dependent on the speed with 
which the individual litmus paper strips respond to the presence of liquid front 
in the sample, especially at times near zero. Thus, the response time is the 
only reasonable explanation that can be offered for the large changes in slope 
shown in the wicking curves for Samples M5, M9, M10, HI, Hll, and H13 with 
water. However, this response time should not affect the final measured value 
of -^m* 
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Figure 62 Freon 113 Wicking Curve for Sintered Nickel Powder Sample Ml 



Figure 63 Freon 113 Wicking Curve for Sintered Nickel Powder Sample M4 
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64 Freon 113 Wicking Curve for Sintered Nickel Fiber Sample H3 
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Figure 65 Freon 113 Wicking Curve for Sintered Nickel Fiber Sample H6 
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The wicking characteristics of the porous materials under investigation are 
summarized in Table 7. The final equilibrium heights for the wicking tests 
utilizing water and Freon 113 are listed. The values listed show that the cap- 
illary forces in the sintered screen samples are relatively weak, while those 
associated with the sintered fiber and powder samples can be quite high, due 
presumably to the finer pore structure obtainable in the latter groups. 


TABLE 7 

Wicking Characteristics of Porous Materials 

(£ m ). water (£ m ) Freon 113 

Sample No. inches inches 


Ml 

9.9 

2.4 

M2 

9.7 

— 

M3 

>15. 9 

— 

M4 

>16.0 

>7. 0 

M5 

8.4 

— 

M6 

6. 9 

2. 0 

M7 

1.9 

— 

M8 

>4.4 

— 

M9 

6. 9 

— 

M10 

9.2 

— 

HI 

>15. 9 

4.2 

H2 

>16. 1 

4. 9 

H3 

>16.1 

8.2 

H4 

>16.0 

— 

H6 

>16. 0 

3.6 

H7 

>16.0 

— 

H8 

>15.2 

— 

H9 

>16.4 

— 

H10 

>16.2 

— 

Hll 

6.1 

— 

H12 

9. 0 

— 

H13 

5.3 

— 


Temperature: approx. 75°F 
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The difference between equilibrium heights for the two liquids can be explained 
in terms of the density, surface tension, and contact angle. As shown earlier, 
the minimum radius of curvature that a wick can support is given by 


^min 


2go<r 


(81) 


where Umin 

= minimum radius of curvature, and it is assumed that 
the meniscus is a segment of the surface of a sphere 

a 

= surface tension 


= density of the liquid 

go 

= proportionality constant in Newton's second law 

g 

= local acceleration of gravity 


The radius of the pore in which the minimum radius occurs is 


r p l*min cos ^ ( 82 ) 

where <p is the contact angle between the liquid and the solid material of the 
wicking sample. 


If it is assumed that the minimum radius of curvature occurs in the same size 
pores in both the water wicking tests and the Freon 113 wicking tests, then 
Equations (81) and (82) can be combined to arrive at an expression relating the 
equilibrium heights, given by 


(-^m)w _ / °~ w \ / ( ^L)F \ / (cos )w \ 

(^m)F \ a F / \ ( / \ ( cos ^F / 


(83) 


where the subscripts w and F refer to water and Freon 113 respectively. 

If the contact angles associated with water and Freon 113 are assumed to be 
equal, and values of fluid properties are substituted, then Equation (83) 
becomes 




5. 94 


(84) 


However, when the experimental values tabulated in Table 7 are used to calcu- 
late the equilibrium height ratio, the resulting values listed in Table 8 are less 
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than those given by Equation (84). This would indicate that the contact angles 
in the water tests and the Freon 113 tests are not equal. 

TABLE 8 


Comparison of Water and Freon 113 Wicking Tests 


Sample No. 


( % m)w 

( 2 ) 

\ m ; F_ 


meas- 

ured 


0 W ( calculated) , 

0-p, - 0 Eq. 83 
r degrees 


(^m)w ( calculated) , 
Eq. (84) 
inches 


Ml 

4.1 

46 

14.3 

M6 

3.4 

55 

12.0 

HI 

>3.8 

< 50 

24. 9 

H2 

>3.3 

< 57 

29.2 

H3 

>1.96 

< 70.8 

48. 7 

H6 

>4.4 

<42 

21.4 


Freon 113 itself is an excellent solvent, which means that it has a contact angle 
of zero or near zero in a relatively clean wick. If it is assumed that 0jr = 0°, 
the contact angle of water in the wicking samples can be calculated using 
Equation (83) and the experimental values of the equilibrium height ratio. These 
calculated water contact angles which are listed in Table 8 show that the contact 
angle with water could have been as high as 55 degrees in some of the wicking 
materials. 

Contact angle data given in Figure 38 shows that the contact angle formed be- 
tween water and a clean pure nickel surface is about 42 degrees, while that 
formed between water and a clean nickel oxide surface is zero degrees. If the 
assumption used above that the minimum pore sizes are equal in both the water 
and Freon 113 tests is correct, then the non-zero value of 0 W listed in Table 8 
results from 1) a wick that was not truly clean, or 2) a wick in which the nickel 
was not completely oxidized. As outlined in Section V C , every effort was 
made to keep the samples clean. Also, each sample was air dried several times 
in an oven at 500°F for two hours, and appeared to be well oxidized (evidenced 
by a dark color and a dull finish) after each oven exposure. Thus no answer can 
be given as to which of the reasons, 1) or 2) above, was responsible for the 
apparent non- zero contact angles tabulated in Table 8. 

Based on the Freon 113 data, the water equilibrium height that would have oc- 
curred had both fluid contact angles been zero can be calculated from Equation 
(84). These values are listed in Table 8 and can be seen to be significantly higher 
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than the experimental values given in Table 7. However, to evaluate each 
wicking material for use in a vapor- chamber fin experimental model, only the 
experimentally determined values of Jl m listed in Table 7 should be used. 

This is the most realistic approach to take because the wicking material that is 
used in the fin will be made and treated in the same way that the wicking sam- 
ples were. Under these conditions, the experimentally determined value of 
JL m for a given fluid is the quantity that should be used in any fin design calculations. 

E. Summary 

To obtain data on the minimum liquid- vapor interfacial radius of curvature that 
a wick would support, tests were conducted to determine the final equilibrium 
height to which a liquid would rise in a wicking material. The wicking materials 
were chosen from the porous material samples listed in Section IV. Two liquids, 
water and Freon 113, were used in the wicking rise tests. 

In general it was found that the time necessary to reach an equilibrium height 
can be quite long, i. e. , in the order of months. In some of the tests conducted, 
a true equilibrium height was not reached, and the test was terminated when the 
liquid front velocity fell below 0. 3 inch in 3 days. 

In the tests involving water, it was found that capillary forces in the best wick- 
ing materials caused the water to rise rapidly to the top of the sample (16 inches). 

In poorer wicking samples the water front height was as low as 1. 9 inches. 

Wicking rise tests involving Freon 113 yielded equilibrium heights lower than 
those utilizing water, due primarily to the lower surface tension of Freon 113. 
Comparison of tests involving both liquids, showed that despite the oxidation 
and cleaning procedures followed, the contact angle in the water wicking tests 
could have been as high at 55 degrees. 
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VI. PERMEABILITY TESTS 


A. Objective 

The object of the permeability tests was to determine the wick friction factor 
for several of the porous materials described in Section IV, and, in particu- 
lar, to determine the dependency of Kq on flow rate, temperature, fluid proper- 
ties and time. 

The wick friction factor range of the metal wicks used in this study is compared 
to those of other porous materials in Table 9. 

As defined in Equation (25), K-^ is the reciprocal of permeability, or 


1 permeability 
with the units of (length) -2 . 


TABLE 9 

Flow Resistance of Various Porous Media 


Substance 

cork board 
hair felt 
cigarette 

metal wicks (P&WA data) 
sand (loose beds) 
Fiberglas 
soils 

bituminous concrete 
sandstone (oil sand) 
agar-agar 
leather 

black slate powders 
brick 

silica powder 
limestone, Dolomite 


Wick Friction Factor 
l/ft 2 x 10~ 8 


0.62 to 2. 8* 

0.77 to 1.1* 

0. 84* 

0. 80 to 61 
5. 2 to 47* 

18 to 39* 

66 to 3,200* 

40 to 9,300* 

310 to 1,900,000* 
2,100 to 47,000* 

7, 700 to 9, 800 * 
7,700 to 19,000* 
4,200 to 190,000* 
18, 000 to 72,000* 
21, 000 to 470, 000 * 


* Scheidegger, A.E. , The Physics of Flow through Porous Media, Macmillan, 
1960, p 87 
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Ki can be measured by passing fluid at a known rate of flow through a rectangu- 
lar wicking sample of known size. The pressure drop across the sample can 
be measured, and can then be calculated from the expression 

where A T is the total cross-sectional area of the sample perpendicular to the 
liquid flow direction, Ax is the length of the sample, PL and pL are the liquid 
density and viscosity, AP-f is the pressure drop across the length of the 
sample, and Wl is the liquid flow rate. 

As mentioned in Section III D , K, is dependent on the geometry and dimensions 
of the passages in the wicking material. 

At the low rates of liquid flow encountered in the vapor- chamber fin, Darcy's 
law? for laminar flow through porous media will apply, meaning that should 
be independent of the flow rate and liquid properties and therefore of liquid 
temperature. 

The permeability investigations of Ginwala, Blatt, and Bilger^ show that in their 
work is initially constant with time, but increases in magnitude after a cer- 
tain amount of flow has passed through the wicking material. This phenomenon 
seems best explained by assuming that tiny gas bubbles are present in the liquid 
in the wick. If these bubbles stick in the matrix of the wicking material, they 
might decrease the effective flow area in the wick and increase K-^. The same 
thing could have been true if foreign particles were present in the liquid. The 
permeability apparatus described in the next section was designed to eliminate 
these problems. 

B. Description of Permeability Apparatus and Test Specimens 

Each test specimen of wicking material was 2 inches wide and 3 inches long in 
the flow direction. The porous material was approximately 0. 1 inch in height 
and was bonded to a backup plate 0. 025 inch thick made of the same material as 
the wick sample. Figure 66 shows a cross-section of a wick sample normal to 
the flow direction. The flow entrance and exit sections were cut by the electro- 
chemical machining process (described in Section IV C ) to prevent any wick 
material distortion. 

A schematic drawing of the permeability apparatus is shown in Figure 67. The 
apparatus consisted of a closed loop fluid system with a mechanical pump which forced 
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Figure 66 Dimensions of Permeability Samples 
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Figure 67 Schematic Diagram of Permeability Apparatus 




a liquid through the wick. The wick was enclosed by gaskets in a sample holder 
so that the fluid flowed through the wick. 

The tubing was 1/4 inch 316 stainless steel tubing, chosen for its corrosion-resis- 
tant properties. The valves used throughout the system were made of 316 stain- 
less steel also. The filter was made of 316 stainless steel sintered powder metal 
and had a five-micron filter rating. The filter prevented microscopic foreign 
material from lodging in the wick test sample. The heaters and coolers used 
were constructed of stainless steel material to prevent corrosion in the system. 

The accumulator tank was a standard 2 l/2-gallon high pressure hydropneumatic 
transfer barrier accumulator. During operation, the working fluid was only in 
contact with the Buta-N bladder or the stainless steel inserts and fittings asso- 
ciated with the bladder. Any desired pressure level could be set in the system, 
and yet a separation could be maintained between the pressurizing gas and the 
working fluid. This separation of the gas and liquid prevented the gas from 
dissolving in the liquid, thus maintaining a gas-free system. 

The pump for the system was a sealless positive-displacement gear pump. This 
type of pump prevented leakage of air into the system and working fluid leakage 
from the system. The impeller gears were made of Teflon and the housing 
material of stainless steel. The drive between motor and pump was coupled 
magnetically through a solid wall. 

Photographs and a sketch of the wick sample holder are shown in Figures 68, 

69, and 70. The holder was constructed of 316 stainless steel and contained 
the necessary fittings for measuring the temperature of the working fluid and 
the pressure drop of the fluid through the wick. The system was designed to 
allow for wick holder removal while preventing the leakage of air into the sys- 
tem. Upon installation of the wick holder, the holder assembly could be evacuat- 
ed to minimize the transfer of entrapped air in the holder to the working fluid. 
Gasket material of neoprene was used to prevent leakage of the fluid around the 
wick sample. 

The manometers chosen for measuring the pressure drop through the wick sample 
were of the thirty- inch U~tube variety. Readings to the nearest 0. 1 inch and 
estimates to the nearest 0. 05 inch of fluid differential could be made. These 
manometers contained air bleeds in each leg for removal of entrapped gas in the 
lines. Manometer taps (see Figure 70) were located on the wick holder to obtain 
a pressure drop measurement within the wick, thus removing any entrance or 
exit effects. Two other manometers measured the pressure drop across the 
entire wick sample. 


PAGE NO. 97 



° N 33V d 






Figure 6|9 Sample Holder for Permeability Apparatus 


RASE NO. 99 


2 IN. MANOMETER TAP 


2 IN. MANOMETER TAP 



Figure 70 Sketch of Wick Permeability Sample Holder Showing Pressure Taps 

Each manometer fluid was chosen both for its specific gravity and its immisci- 
bility with the test liquid. When water was the test liquid, the manometer fluid 
used was acetylene tetrabromide (specific gravity 2. 962 at 70°F) to measure the 
higher pressure drops, and diphosgene (specific gravity 1. 6175 at 79°F) to 
measure the lower pressure drops. When Freon 113 (specific gravity 1. 572 at 
75°F) was used as the working fluid, two different solutions of thallium malonate- 
thallium formate (specific gravities 2. 0151 and 3. 6437 at 75°F) were used to 
cover the pressure-drop range. 

The flowmeters chosen for flow measurement were Fisher-Porter Tri-Flat 
variable-area meters. Three meters were selected to cover the flow ranges. 

The ranges for these meters were 0. 2 to 1. 4 lbm/hr, 0. 5 to 15 lb m /hr and 

1. 5 to 18 lb m /hr for water; and 0. 3 to 1. 8 lbm/hr, 1. 0 to 18 lbm/hr, and 

2. 0 to 25 lb m /hr for Freon 113. 

Temperature measurements were made with chromel-alumel thermocouples con- 
nected to a Honeywell Brown Electronik readout instrument. The accuracy of 
this system is estimated to be ±2°F. Pressure measurements were made using 
Bourdon- tube type gages with a maximum scale reading of 30 psig. Scales were 
in 0. 2 psi increments and readings could be estimated to the nearest 0. 1 psi. 
These gages had bleed ports for removal of entrapped gases in the lines. Experi- 
mental pressures were in the range of 3 to 15 psig. 
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Also included in the permeability apparatus was a means of ensuring that the 
working liquid was not saturated with gas. In preliminary tests, deaerating 
the liquid was done in the closed loop of the permeability apparatus by raising 
the liquid temperature with a resistance heater, causing the solubility of the 
liquid to decrease. The gas expelled from the liquid was collected and bled 
from a high point of the apparatus. The liquid was then cooled to ambient 
temperature, and was in a partially deaerated state. 

Previous experimenters found that increased with time for a fixed flow rate 
and theorized that air, previously absorbed in the liquid, became trapped in the 
wick as bubbles. The trapped air blocked some flow passages and increased both 
the pressure drop and the apparent friction factor K-^. With deaerated liquid, 
the tendency to form bubbles should be reduced. In order to determine the capa- 
bilities of the apparatus to provide values of independent of time, the pre- 
liminary tests were performed with water flowing through a stainless steel sintered 
fibrous metal wick. These tests indicated that during a run of 15-1/2 hours at 
a constant flow rate of 12. 9 lb m /hr, Kq increased 25 per cent. Since other in- 
vestigators ^found that for water saturated with gas the flow through wicks was 
practically stopped during an endurance run, the deaerating method appeared to 
be an effective solution of the problem. However, this method of deaerating the 
liquid in the closed loop is limited in effectiveness because all of the fluid in the 
apparatus cannot be deaerated. The fluid in the pressure gage lines, U-tube lines, 
and the accumulator tank (see Figure 67) is not deaerated because it does not 
flow in the closed loop. 

In an effort to make K-^ constant with time, an additional deaerating apparatus 
was made. This apparatus was separate from the permeability apparatus, and 
its purpose was to supply completely deaerated liquid for filling the permeability 
apparatus. This new apparatus was simply an air-tight vessel made from stain- 
less steel and equipped to heat and boil about five gallons of liquid. Figure 71 
shows the vessel in front of the permeability apparatus during filling. Although 
the filling procedure is explained in more detail below, the basic approach was 
as follows. The test liquid was put into the vessel and brought to a boil. The 
gas was expelled from the liquid due to the high temperature and bled from the 
vessel. This deaerated liquid was then passed into the evacuated permeability 
system. The degree of gas saturation of the water in the permeability apparatus 

after filling with liquid from this vessel was estimated to be less than 20 per 

* 

cent. 

Liquid treated in this way yielded no gas when it was passed through the deaerat- 
ing section of the permeability system loop. Hence, a definite improvement in 


PAGE NO. 101 



decreasing the degree of gas saturation was made by the addition of the deaerat- 
ing vessel. After filling the system with water from the new vessel, an endurance 
test run with the same sample and flow rate resulted in no increase of for a 
period of 25 hours. Thus the time- dependency of K-^ found by other investigators 
appears to have been due to the presence of dissolved gas in the test liquid. 

C. Test Procedure for Permeability Apparatus 

1. Filling of Permeability Apparatus 

The following procedure describes the filling of the permeability apparatus 
using the deaerating vessel shown in Figure 71. Because the permeability 
accumulator tank holds about 2-1/2 gallons of liquid in reserve, this filling 
procedure was only necessary after five or six samples were tested. 

1) The permeability apparatus was drained of any fluid. This involved 
taking the apparatus apart at several locations and draining low 
points. 

2) Nitrogen gas was passed through the apparatus to evaporate any 
remaining fluid. 

3) A vacuum pump was attached to the apparatus so that the entire 
system (without the wick holder) and the tubing connecting the 
apparatus to a closed liquid outlet valve on the deaerating vessel 
were evacuated to a pressure no greater than 1 mm of Hg. 

4) The deaerating vessel was filled with the test liquid, and the liquid 
was then brought to its saturation temperature. During this heat- 
ing process the valve at the top of the vessel was opened for expell- 
ing both gas, which was expelled from the test liquid, and test 
liquid, which drained out because of its expansion. Eventually 
only test fluid vapor came from the top valve on the vessel and a 
pressure could be built up in the vessel by closing the valve. A 
pressure between 2 and 5 psig was maintained by turning the heater 
on and off as a control. 

5) With a pressure less than 1 mm of Hg in the permeability apparatus 
and a pressure of about 5 psig in the vessel, the liquid outlet 
valve on the vessel was opened to begin filling of the permeability 
system. During the process of filling, the evacuation was continued 
at the highest point of the apparatus. When the test fluid in liquid 
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form appeared at the high point, the vacuum pump was disconnected 
from the apparatus. A half hour or more was required after initially 
opening the liquid valve to fill the apparatus. After the pressure 
of the apparatus became the same as the pressure in the vessel, 
the apparatus was valved off from the vessel. To fill the system 
with liquid usually required two-thirds of the contents of the vessel. 

2. Wick Permeability Testing Procedure 

The wick permeability testing procedure was as follows: 

1) The permeability apparatus (except for the sample holder) was 
filled with deaerated liquid as described above. 

2) Measurements of the wick dimensions necessary for evaluating K ^ 
were made. 

3) The wick was cleaned using the same procedure as for the wicking 
rise samples (Section V c ). 

4) The wick and sample holder were assembled as shown in Figure 70. 

5) The sample holder was filled with deaerated liquid from the rest of 
the permeability apparatus by evacuating the holder to at least 

1 mm of Hg and then opening the inlet and outlet sample holder 
valves. The evacuation was discontinued when the holder was full 
of liquid. 

6) The U-tubes were bled, using the positive pressure in the perme- 
ability apparatus. 

7) The test program performed with the permeability apparatus in- 
cluded the evaluation of K-^ for 13 different wicks, evaluation of 

at various flow rates for 12 wicks, the evaluation of at 

various temperatures for three wicks, the evaluation of for 

long testing periods for two wicks, and the evaluation of Kj using 
a different test fluid for four samples. 

D. Test Results 

The permeability results obtained with water on thirteen samples selected 
from those listed in Section IV A are presented in Table 10 and Figures 72 
to 78 in terms of the wick friction factor K^. 
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The porosity values listed in Table 10 for the samples tested may be slightly 
smaller than the porosity values given in Section IV for the respective samples. 
This is the result of a slight compressing of the softer permeability samples 
when being installed between the gaskets of the wick holder. A slight deflection 
of the gaskets on the top and bottom of the sample was necessary to prevent any 
leaks, and on the softer porous materials the pressure associated with the 
gasket deflection caused some compressing. The effect of compression has 
been accounted for in the porosity and permeability results presented in this 
section, by using the sample dimensions measured after each permeability test. 

TABLE 10 

Wick Friction Factor for Porous Material Samples 
Test Fluid: Water 

Wick Friction Factor 

Sample No. Porosity, % x 10~8, l/ft^ 


HI 

83.6 

21. 2 

H2 

81.4 

27.6 

H3 

68.9 

61.3 

H6 

86. 9 

30.2 

Hll 

88.7 

1.70 

H12 

80.7 

4.75 

H13 

81. 8 

0. 80 

M2 

64. 5 

3.4 

M4 

54. 0 

11.5 

M6 

69.6 

3.1 

M7 

62.5 

1.4 

M8 

58.5* 

6.1 

M10 1 

65. 1 

12. 0 


* Significant amount of change in porosity due to compression 

In the range of flow rates encountered in the vapor- chamber fin, Darcy’s law 7 
for flow through porous media applies, meaning that as given by Equation 
(86) should be independent of flow rate. Figures 72 to 75, showKj vs flow rate 
per unit area (i. e. , the flow rate w^ divided by the cross-sectional area At) 
for different sintered screen samples, sintered powder samples, and sintered 
fiber samples. These figures demonstrate that for a certain range of flow rate 
per unit area, is independent of flow rate. 
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Figure 72 Wick Friction Factor vs Flow Rate per Unit Area for Sintered 
Nickel Screen Samples 
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Figure 74 Wick Friction Factor vs Flow Rate per Unit Area for Sintered 
Nickel Fiber Samples 
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FLOW RATE PER UNIT AREA, W/A T x|0‘ 3 

LBm /HR FT 2 

Figure 75 Wick Friction Factor vs Flow Rate per Unit Area for Sintered 
430 Stainless Steel Fiber Samples 
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In interpreting these figures it should be noted that the points given at very low 
rates of flow per unit area become inaccurate for either or both of two reasons. 
First, the instrument for measuring flow rate for flow rates per unit area less 
than 1 x 10 3 lb m /hr ft 2 (experimental flow of 1. 4 ijigl ) becomes increasingly 
inaccurate as flow rate decreases. Second, the inaccuracy in measuring the 
pressure drop for samples with relatively low values of Ki becomes greater 
at low flows because of the small pressure drops. These two reasons account 
for the scatter of points at the lower flow rates per unit area. 


The data presented in Figure 72 for three sintered screen samples, demon- 
strates the independency of Ki on flow rate over most of the range of flow rates 
per unit area tested. In Figure 73, in which data for sintered powder samples 
is presented, the value of K for sample M4 remains constant from 1 x 10 3 to 
9 x 10 3 lb m /hr ft 2 , but then increases slightly from 9 x 10 3 to 15 x 10 3 
lb m /hr ft2. This increasing K-^ indicates that the Darcy flow regime has been 
exceeded. Within the Darcy flow regime friction is the only cause of pressure 
drop, but at higher flow rates inertia effects resulting in drag forces start to 
enter as a pressure drop component, causing to increase. This upward 
trend of K-^ at flow rates per unit area above 9 x 10 3 lb m /hr ft 2 is also present 
in Sample M10 of Figure 72, Sample M2 of Figure 73, and Sample H12 of 
Figure 75. All of the samples tested over the flow rates per unit area of about 
1 x 10 3 to 9 x 10 3 lb m /hr ft 2 resulted in values of having very little or no 
dependency on flow rate. 


Figure 76 is a plot of wick friction factor vs fluid temperature of several dif- 
ferent samples for tests with water. In the temperature range of 75 to 160°F the 
friction factor remained constant while the liquid viscosity changed by a factor 
greater than two. This shows that, in the range tested, as given in Equation 
(86) is independent of the variation of fluid properties with temperature. 

Figure 77 is a plot of wick friction factor versus time for wick Sample H12 with 
deareated water and M2 with aerated water. Over a test time of 84 hours there 
was no change in the value of for Sample H12 with deaerated water as the 
test fluid. Over a test time of 53 hours, the value of increased 18 per cent 
for Sample M2 with aerated water as the test fluid. This trend was experienced 
by other investigators 11 with test liquids which had not been deaerated. For 
all permeability test data presented previously in this report, deaerated test 
fluid was used. Thus, it is expected that the results presented are constant 
with respect to time. 

Figure 78 is a plot of wick friction factor versus flow rate per unit area for 
four different samples with two different test fluids, water and Freon 113. 
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Figure 76 Wick Friction Factor vs Temperature for Sintered Metal Wicks 
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Figure 77 Permeability Duration Runs Using Aerated and Deaerated Water for 
Two Sintered Metal Wicks 
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Figure 78 Wick Friction Factor vs Flow Hate per Unit Area with Test Fluids 
Water and Freon 113 for Sintered Metal Wicks 


PAGE NO. 113 





In general the variation of the Freon 113 data points from the average value of 
Ki obtained from the water results was within the probable uncertainty discussed 
in the next section. Larger variations than expected occurred for only three 
Freon 113 data points on Sample H-13 in the Darcy flow regime. As Darcy's 
law implies, seems to be independent of the test fluid properties for flow 
within the Darcy flow regime. 

E. Accuracy of Permeability Results 

To determine K-^ as defined in Equation (86), pressure-drop data was taken for 
two lengths of each sample: the entire sample length of about three inches, 
and an internal length of about two inches. The nominal three-inch dimension 
of the sample was measured within ±1 per cent. However, the electrochemically- 
cut entrance and exit sections of the wicks were not always perpendicular to the 
nominal three-inch dimension, and for some samples an angle of 10 degrees 
from the perpendicular existed, caused by the cutting process. Because this 
discrepency might have caused a slight flow maldistribution at the sample en- 
trance and exit, data from both lengths was compared in Figure 79, demon- 
strating a usual variation of less than ±3 per cent. The nominal two-inch dimen- 
sion is fixed by the construction of the pressure taps in the sample wick holder, 
as shown in Figure 70. The actual measurement was 2. 06 inches for the center- 
to-center distance. 

The probable uncertainty of determined with water in the permeability ap- 
paratus is summarized in Table 11. The table includes the expected accuracy 
of all parameters used in calculating from Equation (86). All of the results 
of presented in this report were calculated using pressure drops measured 
from the two inner pressure taps, except where the pressure differential between 
these taps read below two inches of manometer fluid. For this latter case the 
pressure drops between the outer set of pressure taps (3 inches center-to- 
center distance) was used. The data from the inner set of pressure taps was 
preferentially used in comparing the different samples to eliminate end effects. 
From Table 11 the expected accuracy of Ki presented in this report may be 
expressed by dividing the samples into two groups. For samples with K^< 

4. 0 x 10 8 l/ft 2 the probable uncertainty of the Kj results are within 9. 8 per 
cent for flow rates per unit area above 4. 4 x 10^ lb^/hr ft^, and within 5. 3 per 
cent for flow rates per unit area above 7. 3 x 10* 5 lbm/hr ft . For the samples 
with > 4. 0xL0^ l/ft 2 , the probable uncertainty of Ki is within 9. 8 per cent 
for flow rates per unit area above 1. 4 x 10^ lb m /hr ft 2 , and within 5. 3 per 
cent for flow rates per unit area above 4. 1 x 103 lb m /hr ft 2 . All experimental 
data obtained for water except one point was within the probable uncertainty 
range from the mean line through the data. 
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Figure 79 Wick Friction Factor for Sample M4 vs Flow Rate Using an Internal 
and External Sample Length 

F. Summary 

The wick friction factor K of each of thirteen of the samples listed in Section 
IV was measured and was found to range from 0. 80 x 10^ l/ft 2 for a high- 

CS O 

porosity sintered fiber sample to 61.3 x 10° 1/ft for a low-porosity sintered 
fiber sample. 


The wick friction factor values were determined in a range of flow rates where 
Kj was independent of the flow rate, i. e. , in the Darcy flow regime. 

The wick friction factor values were found to be independent of temperature. 


The wick friction factor was found to be independent of time if the test fluid 
was deaerated. However, an increase of with time was found if the test 
fluid was aerated. 


The value of Ki appeared to be independent of the test fluid properties as deter- 
mined from tests on samples using both water and Freon 113. 
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TABLE 11 


Accuracy of Wick Friction Factor, 


Parameter 




X 



or 







Dimension 

Accuracy of 3-Inch Results 

Accuracy of 2- Inch Results 




Per cent 



Per cent 


Magnitude 

Uncertainty 

Uncertainty 

Magnitude 

Uncertainty 

Uncertainty 

density (H 9 0 at 75°F) 

62.2 <±0.1 

<±0.16% 

62.2 

<±0. 1 

<±0.16% 

ibjk 3 



±1.78% 



±1.78% 

viscosity (H 2 O at 75°F) 
lb m /ft sec x 1 0 3 
length, Ax 

0.616 

±0.01 

0.616 

±0.01 



3.00 1 

±0.020 

±0. 67% 

2.060 

<±0. 01 

<±0.48% 

inches 




0.20 1 



cross-section area, A c 

0.20 1 

±0. 002 

±1. 00% 

±0. 002 

±1. 00% 

(inches) 2 







flow rate 

Flowmeter 

Reading 






flowmeter No. 1 







0. 41 lb m /hr 
1. 24 lb™/hr 
flowmeter No. 2 

10.5 

±0. 1 

±1. 68% 2 




18.8 

±0. 1 

±2.25% 2 

Same As Accuracy 





of 3' 

" Results 


1. 24 lb m /hr 

5.0 

±0.25 

±7. 24% 2 




3. 84 lb m /hr 

9.5 

±0.25 

±3. 81% 2 




12.2 lb^ /hr 
m 

22.0 

±0.1 

±1. 15% 2 




pressure drop 

U-Tube 



U-Tube 




Reading 



Reading 



0. 056 psi 

2.50” 

±0. 05" 

±2. 00% 

0.79" 

±0. 05" 

±6. 32% 

0. Ill psi 

5. 00" 

±0. 05" 

±1. 00% 

1.57" 

±0.05" 

±3. 18% 

0. 666 psi 

30.0" 

±0. 05" 

<±0. 50% 
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TABLE 11 (Cont'd) 


Parameter 

or 

Dimension 


Accuracy of 3-Ineh Results 

Percent 

Magnitude Uncertainty Uncertainty 


Accuracy of 2-Inch Results 

Percent 

Magnitude Uncertainty Uncertainty 


±0.05 
±0. 05 

wick friction factor, 

Kl 

for all flows 
from 0. 41 lb m /hr 
to 12. 2 lb m /hr 
and for AP 

of 0. 056 psi to 2.12 psi 

for flows 3. 84 lb m /hr 
and AP 0. Ill psi 


pressure drop 

0. 666 psi 
2. 12 psi 


U-Tube 

Reading 

9.4" 

30. 0" 


±0.53% 9.4" ±0.05" ±0.53% 

±0.50% 30.0" ±0.05" ±0.50% 


±7.75% 3 


±9.75% 3 


±4. 26% 3 


±5. 27% 3 


The exact magnitude varied slightly in different samples 
2 These uncertainties take into account the instrument accuracy and the calibration error 
3 These uncertainties are obtained from the square root of the sum of the uncertainties squared 
for the parameters in Equation (86) 



VII. DISCUSSION OF CAPILLARY PUMPING PARAMETER 


A. The Capillary Pump 

The porous material of a wick acts as a liquid pump because of the action of 
capillary forces. These forces are produced by the interaction of the liquid- 
vapor interface in the wick and the internal solid surface of the wick itself. 
Whenever the capillary forces bring about any movement of the liquid in the 
wick, frictional forces arise that oppose the liquid motion. 

A capillary pump can be evaluated by calculating the maximum flow rate that 
it can achieve. This was essentially done in deriving the analytical model of 
the vapor-chamber fin in Section m. However, to eliminate the complication 
of mass addition along the length of the wick, the following capillary pump 
model will be considered to illustrate what this output depends upon. 

Shown in Figure 80 is a rectangular piece of wicking material acting as a 
capillary pump. The material has a total flow area of At and a flow length 
of Ax. 

At the front face, a wetting liquid is being evaporated at such a rate that the 
wicking material is providing a maximum flow. If this is true, then the pres- 



Figure 80 Model of a Capillary Pump 
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sure drop (AP)^ across the liquid- vapor interface at the plane of evaporation 
is given by the Laplace- Young equation 


<**>r 


2 a 

^min 


(87) 


where <r is the fluid surface tension and R m i n is the minimum radius of curva- 
ture of the liquid- vapor interface that the wick can support. It has been assumed 
that the menisci are segments of the surface of a sphere. The pore radii pro- 
ducing the curvature in Equation (87) will be equal to Rmin onl y if the contact 
angle is zero. 


The pressure drop given in Equation (87) must be balanced by the frictional 
pressure drop caused by the flow of liquid through the wick, so that 

t 


(AP) a - <AP) f (88) 

Combining Equation (88) with Equation (87) and the definition of the wick fric- 
tion factor, Equation (25), yields an expression for (wL)m.ax> the maximum 
flow rate that the capillary pump pictured in Figure 84 can produce , or 


( w l) max 


2 PL_? \ ( Ac V 1 

/ \ Ax /yK^Rmin 


) 


(89) 


where K^» the wick friction factor, is a constant for flow in the Darcy regime. 
Thus, given a wick size and a fluid, the maximum output of the capillary pump 
depends on the parameter 


1 

■^1 ^min 


(90) 


As discussed in Section III D, Ki is dependent only on the pore geometry in the 
wick. However Rmin is a function of the pore geometry in the wick, the wick 
material, and the liquid. The parameter in (90) shows that any evaluation of a 
capillary pump must include both a consideration of the liquid and the wick itself. 


In Section V it was shown that R m i n for a given wick could be determined from 
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R min “ 


_ 2 SQ..f„ 

^L^m 


(91) 


where £ m is the maximum height to which a liquid will rise in a vertical 
sample of the wick when exposed to a gravity field, g. 


By combining Equations (91) and (89) the maximum capillary pump output is 
given by 


where the subscript WR refers to the conditions under which J?m was deter- 
mined. Equation (92) can be used to express the amount of energy "pumped" 
by multiplying both sides by the latent heat of vaporization. 

Hereafter, the term Z m/K ^ will be referred to as the capillary pumping para- 
meter. It should be emphasized that this parameter is a function of the liquid, 
the wick geometry, and the wick material. Equation (92) shows that -£m/K^ 
is a significant parameter for comparing the pumping capabilities of various 
wicking materials of the same size with the same liquid. 

B. Experimental Capillary Pumping Parameter Results 

Values of the capillary pumping parameter for both water and Freon 113 are 
listed in Table 12. Also tabulated are values of porosity, wick friction factor, 
and equilibrium heights taken from the experimental data in Sections IV, V, 
and VI. 

From Table 12 it can be seen that the capillary pumping parameter for Freon 113 
is lower than that of water for each sample tested. 

Although Ki is independent of liquid properties, is not. Values 

ofj m with Freon are lower than those with water for a given sample, for the 

reasons given in Section V D , which results in lower values of^ m /Kq. 

The values of£m/ K l h 1 Table 12 show that sintered fiber samples form the best 
capillary pumps, followed by sintered powder samples. The sintered screen 
samples demonstrated the lowest potential as capillary pumps. 


PAGE NO. 120 


TABLE 12 


Capillary Pumping Parameter 


Sample 

No. 


Porosity, 

% 

<Uv, 

inches 

inches 

K i, 

i/ft 2 

ft 3 

Wf, 

ft 3 

Ml 

powder 

59.7 

9.9 

2.4 



- 




M2 

it 

65. 8 

9.7 

- 

3.40 

x 10 

00 

to 

• 

CO 

oo 

io- 

9 

— 

M3 

M 

47. 7 

>15.9* 

- 


- 

- 



— 

M4 

tf 

54. 0 

>16.0* 

>7.0 

11. 5 

tt 

> 1.16 

it 

>0. 

508 x 10~ 

M5 

tl 

69. 1 

8.4 

- 


- 

- 




M6 

1 1 

69. 6 

6. 90 

2. 0 

3. 10 

it 

1.85 

tt 

0. 

543 " 

M7 

screen 

62. 5 

1.9 

- 

1.40 

tt 

1.14 

it 


- 

M8 

ti 

67. 9 

>4.4 

- 

6. 10 

it ^ 

> 0. 604 

tt 

** 

- 

M9 

tt 

67. 8 

6. 9 

- 


- 

- 



- 

M10 

it 

67. 6 

9.2 

- 

12. 0 

it 

0. 640 

IT 


- 

HI 

fiber 

86. 8 

>15. 9* 

4.2 

21.2 

tt 

> 0.625 

ft 

0. 

165 " 

H2 

tr 

82. 5 

* 

>16. 1 

4.9 

27. 6 

it 

> 0.486 

tt 

0. 

148 " 

H3 

t! 

68. 9 

>16. 1* 

8.2 

61. 3 

tt 

> 0.219 

ft 

0. 

112 " 

H4 

n 

62. 8 

>16. 0* 

- 


- 

- 



- 

H5 

tt 

91. 8 

- 

- 


- 

- 



- 

H6 

tt 

88. 0 

>16.0* 

3. 6 

30.2 

it 

> 0.442 

tt 

0. 

099 11 

H7 

tt 

82. 0 

>16. 0* 

- 


- 

- 



- 

H8 

it 

82. 8 

>15.2* 

- 


- 

- 



- 

H9 

tt 

62. 6 

>16. 4* 

- 


- 

- 



- 

H10 

it 

70. 9 

>16.2* 

- 


- 

- 



- 

Hll 

tt 

91. 6 

6.1 

- 

1. 70 

tt 

2. 99 

ft 


- 

H12 

tt 

80. 8 

9.0 

- 

4. 75 

it 

1.58 

Tl 


- 

H13 

tt 

82.2 

5.3 

- 

0. 80 

tt 

5.54 

tf 


- 

* Indicates that 

water wicked to top of sample used in 

wicking- 

rise 

test 


** Values may have been significantly affected by the sample being compressed 
during permeability test (see Table 10) 
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Also shown by the tabulation of the capillary pumping parameter in Table 12, is 
that Sample H13, having the largest value of£ m /Ki, would serve as the best 
capillary pump, followed by Hll and M2. The table indicates that a good capil- 
lary pump cannot be chosen on the basis of wicking ability alone, for considera- 
tion of frictional forces is equally important. For example, the best material, 
H13, has one of the lowest values of for water, but it also has by far the 
lowest value of Kj. 

The values of porosity listed in Table 12 are the same as those listed in Table 
4 and not those listed in Table 10. The difference between the two tables of 
porosity is due to the fact that some samples were compressed while being 
installed in the permeability rig, thus resulting in the slightly lower values 
listed in Table 10. 

C. Effects of Gravity on Capillary Pumping 

It should be noted that, although H m / K l is a measure of capillary pump output, 
it does not include the effect of gravity or head. If an application involves a 
condition in which a capillary pump must work against gravity, consideration 
must be given to both the value ofJdm/Kl and^ m itself. For instance if a capillary 
pump material is to be used in a situation where liquid must be transported by 
the material over a vertical distance of 10 inches, a material must be chosen 
that has a value of jZ m greater than 10 inches, as well as a high value ofj£ m /Ki. 

This separate consideration of J m in a gravity situation cannot be neglected in 
the case of the vapor-chamber fin. In the analysis of the one-dimensional model 
of the vapor chamber fin in Section III, it was shown that the maximum total rate 
of heat rejection can be obtained from the following equation 

«-■ (-£-) (i rH-fc) « 

where gravity is neglected. 

If Equation (93) is modified to include the effect on the liquid in the wick as it 
flows against gravity, Q max is given by 

W - (^) (-£■ (4r) - 

2 )( ib ) ( If ) (it) <* 4 » 

where 6 is the angle at which the fin is inclined. Equation (94) shows that, for 
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a given fluid and fin size, to maximize the value of Qmax for sin 6 > 0, a wicking 
material must be chosen on the basis of a large value ofjl m /Kj and also a large 
value of Jl m . However, large values of both of these qualities do not occur 
simultaneously in any of the samples studied. As the experimentally determined 
values in Table 12 show, a high value of j£ m /Ki seems to imply a low value of 
■-^m* This might be expected since a porous material that has a large pore size 
will be characterized by low capillary forces (small value ofJ&m), but a large 
pore size is indicative of low frictional forces for a given flow rate (small value 
of K^), all of which results in a large value of £ m /Ki but a low value of 
Section VII D contains a more detailed discussion of the effect of pore size on 
the capillary pumping parameter. 

D. Capillary Pump Pore Sizes 

Some idea of the structure of the porous materials can be obtained by comparing 
the pore size distribution, wick rise and permeability data. 

The most direct way to get the pore size or sizes of a given capillary material 
would be by direct measurement. Such a process is obviously not feasible 
and would result in data difficult to interpret, so that some indirect method of 
measuring the pore sizes in a material must be used. 

The mercury intrusion method of measuring pore size was discussed in Section 
IV. Pore diameters were shown to be given by the equation 


D 


4 <r cos 0 
P 


(95) 


where 

D = diameter of pore 

P = pressure required to force mercury into a pore of diameter D 
cr = surface tension of mercury 

<i> = contact angle between mercury and the wick matrix material 

The median pore diameter, was then defined as the pore size at which 

50 per cent of the connected pore space is filled with mercury. Values of 
D med ta ken from experimental data of Section IV are listed in Table 13. 
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An equivalent pore diameter, D^, can be calculated from the wick friction data 
presented in Section VI. For a given flow rate of liquid through the wicking ma- 
terial, the resulting pressure drop AP over a given sample length Ax can be 
expressed in terms of the Poiseuille equation for flow in a tube, to get 


AP 

Ax 


32 ML u L 



(96) 


where Ml * s liquid viscosity, ul is the average liquid flow velocity through 
the pores of the material and D& is the equivalent diameter defined by this equa- 
tion. As defined here, Dj^ is a variable that contains the effects of tortuosity, 
intereonnectiveness of the pores, and an average pore shape factor and size. 


The above expression can be combined with Equation (86) and the continuity 
equation to obtain 



32 \ 1/2 

f ar k i ' 


(97) 


where is the wick friction factor, defined by Equation (25). The free flow 
area ratio, F^r, is given by Equation (72). In effect, Equation (97) represents 
the porous material as a bundle of straight parallel tubes of equal size. 


The values of Dg listed in Table 13 are based on values of Kq from the mean 
line through data taken using water as the test fluid. These values show that 
the actual numerical agreement between Dr and Dmed is reasonably good at 
least for the fibrous wicks. This is somewhat surprising, for the two methods 
are physically quite different, one using frictional data and the other surface- 
effect data. While each method should be consistent within itself, it would not 
be expected that values of pore size calculated by the two methods would agree 
so closely. This correspondence may mean that in the flow region in which 
was determined, the effects of tortuosity, intereonnectiveness of the pores, and 
the pore shape factor combine in such a way that the assumption of straight 
parallet equal-sized tubes made for Equation (97) is an acceptable simplification. 


Also listed in Table 13 is the equivalent capillary pore diameter computed from 
the wicking data of Section V, using the equation for capillary rise in a tube of 
diameter with a contact angle of zero degrees 


D 


C 


4 g 0 ff 

Pjg£ 

L m 


(98) 
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TABLE 13 


Effective Pore Sizes 


Measured 

Sample Median Pore Equivalent Frictional Equivalent Capillary 


No. 

Type 

Diameter, D me( j 

Pore Diameter, 

OK Pore Diameter, Dq 

Ml 

Powder 

microns 

microns 

Water 

microns 

117 

Freon 

microns 

81 

M2 

Powder 

>70 

116 

120 

- 

M3 

Powder 

- 

- 

< 73 

- 

M4 

Powder 

0 

00 

1 

69.2 

< 73 

<28 

M5 

Powder 

- 

- 

138 

- 

M6 

Powder 

>80 

117 

168 

96 

M7 

Screen 

- 

184 

609 

- 

M8 

Screen 

>100 

91. 3* 

< 263 

- 

M9 

Screen 

- 

- 

169 

- 

M10 

Screen 

— ' 50 

61. 7 

126 

- 

HI 

Fiber 

43. 8 

41.0 

< 73 

46 

H2 

Fiber 

35. 0 

36.4 

< 72 

40 

H3 

Fiber 

26.5 

26. 5 

< 72 

23 

H4 

Fiber 

- 

- 

< 73 

- 

H5 

Fiber 

- 

- 

- 

- 

H6 

Fiber 

40.7 

33. 6 

< 73 

54 

H7 

Fiber 

- 

- 

< 73 

- 

H8 

Fiber 

- 

- 

< 77 

- 

H9 

Fiber 

- 

- 

< 71 

- 

H10 

Fiber 

- 

- 

< 72 

- 

Hll 

Fiber 

> 80 

140 

190 

- 

H12 

Fiber 

>70 

88. 1 

129 

- 

H13 

Fiber 

>80 

213 

219 

- 


* Value may have been significantly affected by the sample being compressed 
during permeability test (see Table 10) 


where (L m is the equilibrium height to which the liquid (water or Freon 113) 
will rise in a vertical sample. 

The discrepancy between the values of Dq for water and Freon 113 for a given 
sample can be explained by assuming that the contact angles for each liquid 
were different. Table 8 gives estimates of the water contact angle assuming 
the contact angle for all Freon 113 tests was equal to zero degrees. 
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The agreement between the values Dq and the values of Dnied and Dj^ * n Table 
13 is reasonably good except in the case of the screen Samples, M7, M8, M10, 
powder Sample M4, and fiber Samples H6 and H12. This is curious since it 
would be expected that for a given material , D c would be smaller than D me( j. 
The latter is a measure of the most populous pore size , while Dq is calculated 
from the equilibrium height, J2 m , which should be a measure of some smaller 
pore size in the wicking material. 

The agreement between these three diameters is particularly good for the fiber 
samples. It appears from this agreement for these samples that the actual liq- 
uid levels recorded in the wicking tests were characterized by a pore size close 
to Dmed- If the smaller pores in the wick are considered to be interconnected 
only by median-size pores, the value of Stm would be expected to depend on the 
median pore size rather than some smaller size since, during a wick rise test, 
liquid would be unable to bridge the median pores to reach smaller pores. Such 
a situation would result in equilibrium heights characterized by D me( j* Figure 
81 is a plot showing the percent variation of Dq from Dj^ for the wicks tested. 
This plot shows that, except for the screen samples, Dq was within ± 60 per 
cent of D-£ for both the water and Freon 113 results. Also, Dj^ was found to be 
close to Dmed f° r the two groups of powder and fiber samples. The large dif- 
ference between Dq and D K for the screen samples might be caused by one or 
both of two factors: 1) may be much smaller than Dq because of a large 
tortuosity for liquid flow in the wick, or 2) a regular occurrence of large pores 
in all of the flow passages of the screens due to the stacking of the screen layers 
might prevent the liquid from rising in a wicking rise test as high as it would 
in a more homogeneous porous material. The resulting value of D c may be 
larger than D K . The smaller discrepancies of Dq and Dj^ for fiber and powder 
samples might be caused by three similar factors: 1) Dk may be smaller than 
Dq because of the tortuosity for liquid flow in the wick, 2) D^ may be smaller 
than Dj£ because smaller pores may be more important in wick rise tests, or 
3) Dq may be larger than Dj^ because of a nonhomogeneous distribution of 
large pores preventing the liquid in the wicking rise tests from rising as high 
as it would in an homogeneous material. 

Figure 82 shows a plot of vs D K for samples tested in water with 

porosity as a parameter. This plot is intended to show that, in a very general 
way, J?in/Ki seems to increase with pore size for the samples tested. 

An explanation of the trends shown in Figure 82 can be given if the assumption 
is made that Dq is equal to Dj£. This assumption is based on the agreement 
between values of these two diameters shown in Table 13 for the powder and 
fiber samples. Using this assumption, Equations (97) and (98) can be combined 
to yield an expression for as follows 
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Figure 81 Comparison of Diameters Obtained from Permeability and Wick Rise Tests 
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But as shown in Section IV, the free-flow area ratio can be expressed as 


TT =6 

AR 


( 100 ) 


if it is assumed that the wicking material is isotropic. Using Equation (100), 
Equation (99) becomes 


£ m _ g o a 6 

Ki ” 8 g P L 


( 101 ) 


In effect Equation (101) reduces a wick to an equivalent capillary tube. For 
such a situation the capillary forces are inversely proportional to the diameter 
of the tube, (i. e. >(? m ~l/Dj£) whereas the frictional forces are inversely 
proportional to the diameter squared, (i. e. Ki ~ l/Dj^ which leads to the 
linear dependence of .£ m /Kq on the diameter shown in Equation (101). 


Equation (101) is plotted in Figure 82 for values of e equal to 60 and 90 per cent. 
The data is in fair agreement with these plots. The points representing the 
sintered-screen samples lie the furthest from these two curves as would be 
expected from the previous discussion of Figure 81. 


Thus Figure 82 shows that the capillary pumping parameter increases with 
increasing equivalent pore diameter. However, the plot shown in Figure 82 
includes the parameter for frictional and capillary forces, and not gravity forces. 
As pointed out in Section VII C , the evaluation of a wicking material for use 
in a gravity situation must include a consideration of the 2 m of the material 
as well as J2 m /K|. 

E. Summary 

The maximum output of a capillary pump was found to be dependent on the 
capillary pumping parameter Jt m / when zero net body forces are considered. 
Experimentally determined values of the capillary pumping parameter are pre- 
sented for thirteen samples. In applications involving a condition in which a 
capillary pump works against gravity, consideration must be given to both the 
value of the capillary pumping parameter and^ m , the equilibrium height of 
liquid in the vertical wick. 

Effective pore diameters were computed from mercury porosimeter data, 
permeability data, and wick rise data. The capillary pumping parameter was 
found to increase almost linearly with the median pore diameter of the wick. 
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VIE. BOILING STUDIES 


A. Boiling in Wicks 

The limitation on fin operation due to boiling occurring in the evaporator wick 
was discussed in Section HI B. A separate test program was performed to 
study the process of evaporation and boiling of a liquid in a wick. 

A number of investigators have studied boiling at wick covered surfaces, but 
none of their researches are directly applicable to the heat transfer process 
in the evaporator section of the vapor-chamber fin. Allingham and McEntire- 1 - 2 
obtained boiling coefficients for a tube embedded in water-saturated ceramic 
fiber wicking material. However, because their data was taken at relatively 
low temperatures (120-160°F) and low heat fluxes (1000 - 10, 000 Btu/hr ft 2 ), 
no conclusion can be drawn about peak heat fluxes. Critical heat flux data 
for a wick covered surface was obtained by Costello and Redeker 13 , but because 
of the path taken by the vapor generated at their heater surface, little correlation 
can be made between their experimental model and the evaporator of the vapor- 
chamber fin . 

The later work of Costello and Frea 14 is perhaps the most germane to the problem 
of boiling in the fin. By wrapping a cylindrical electric heater in wicking ma- 
terial they were able to obtain pool-boiling critical heat fluxes as high as 1.2 
x 10 6 Btu/hr ft 2 in water, with a 50°F difference between the heater temperature 
and the saturation temperature of the liquid. The data of Costello and Frea 
is both gravity-dependent and geometry-dependent, because the vapor produced 
was not vented through the wicking material, but rose through an opening in the 
wick above the heater. Thus, the question of the ability of the wick to support 
the flow of vapor from and liquid to a region of the heated surface was not conclusively 
answered. However, the extremely high value of the critical heat flux obtained 
in their studies was very encouraging. 

B. Objective 

The objective of the tests was to determine the following boiling heat transfer 
characteristics while boiling from a wick-covered planar surface: 

1) The variation of heat flux, q, with AT gat , the temperature at the 

liquid- wick-solid interface minus the saturation temperature , and 
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2) The maximum heat flux that can be sustained before film boiling 
occurs. 

The two objectives above were to be determined as functions of wick mater- 
ial, boiling fluid (water or Freon 113) , wick orientation with respect to 
gravity, degree of wick submergence, and surface modification. 

Three classes of wicking materials were studied: sintered metal screens, sintered 
metal powders, and sintered metal fibers. The results are interpreted to deter- 
mine limitations on boiling from a wick-covered planar surface. The following 
sections contain a description of the apparatus used to accomplish these objectives 
and the results obtained from the test program. 

C . Description of Boiling Apparatus 

1. Overall Design 

Figure 83 shows a diagram and Figure 84 a photograph of the boiling test 
setup, which includes the primary tank containing the test sample, the 
sample heater, and the guard heater; the auxiliary tank with a heater, and 
a liquid heater located between the primary and auxiliary tanks. A sketch 
of the stainless steel^primary tank is shown in Figure 85, and a photograph 
of the interior is shown in Figure 86. The flat test sample to be evaluated 
was bonded to a resistance heater mounted on copper busbars. The top side 
of the heater was coated with an electrical insulator (aluminum oxide). The 
bottom of this heater was firmly pressed against insulating blocks by means 
of busbar adjusting screws. The insulating blocks, two of which had a Ni- 
chrome guard heater sandwiched between, rested on a Textolite platform 
which was mounted in the stainless steel primary tank. A Plexiglas cover 
and Neoprene gasket were used to prevent spillage during the tests in which 
the primary tank was tilted. The overall dimensions of the inside of the 
primary tank were 10 3/4 x 2 1/4 x 5 3/4 inches. 

The auxiliary tank was 3. 5 by 2 by 0. 5 ft and was made of stainless steel. 

The purpose of this secondary tank was to provide a large liquid-vapor 
surface, so that the change in liquid height with time, as the primary tank 
fluid boiled away, was limited to a negligible amount. In this way, the 
degree of wick submergence was kept constant during a test. A heater was 
provided in the secondary tank to preheat the liquid before entering the pri- 
mary tank. 

A fixture was built for tilting thfe primary tank to determine the effect of 
wick orientation with respect to gravity on heat transfer coefficients. The 
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Figure 83 Schematic Diagram of Boiling Test Setup 
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fixture consisted of two steel plates hinged together at one end and connected 
by supporting rods on the other end. The tank, which was bolted to one of 
the plates, could be tilted up to 90 degrees by varying the length of the rods. 
Figure 87 shows the tank operating at a 45-degree angle. 

2. Sample-Heater Assembly 

The wick boiling samples consisted of a piece of porous nickel or stainless 
steel 2 7/8 by 5/8 by approximately 0. 1 inch thick bonded to a backup plate 
of the same material. The exact thickness of each wick used in the boiling 
tests is given in the table below. The samples were cut from the wicks 
used in the wick rise tests by electrochemical machining. Since a wick 
boiling sample was small in size, there should have been little resistance 
to liquid flow in the wick. Hence, for tests in which the sample was hori- 
zontal with respect to the surface of the water, failure should have been 
due to the onset of film boiling rather than the wicking inability of the sample. 

Wick Thickness of Samples Used in Boiling Tests 

Mean 

Wick Sample Thickness, Inches 


M2 - 

sintered nickel powder 

0. Ill 

M7 - 

sintered 50-mesh nickel screen 

0. 106 

M9 - 

sintered 150-mesh nickel screen 

0.080 

H6 - 

sintered nickel fiber 

0.101 

Hll - 

sintered stainless steel fiber 

0. 098 

H13- 

sintered stainless steel fiber 

0.100 


The heater was made of Inconel X for the nickel samples, and stainless 
steel for the stainless samples. This material selection reduced differential 
expansion between the heater and sample which could break the sample- 
heater bond. The top side of the heater was plasma arc-sprayed with a 
0. 002-inch coating of aluminum oxide. A thin layer of epoxy bonded the 
coated surface to the sample backup plate. A 1/8 inch Transite block was 
placed under the heater and held in place by two l/8-inch wide layers of 
epoxy along the side edges of the heater and blocks. Figure 88, which is a 
photograph of a completed sample, shows these edges of epoxy; and Figure 
85 shows a schematic of sample-heater construction. This type of sample- 
heater construction had the following advantages ; 

a. The aluminum oxide coating provided the electrical resistance 
necessary to prevent current leakage across the sample, while 
the thermal conductivity (10 Btu/hr ft°F) remained quite high for 
an electrical insulator. Hence, the temperature drop across the 
aluminum oxide was negligible for the heat fluxes of interest 
(10, 000 Btu/hr ft 2 to 150, 000 Btu/hr ft 2 ). 
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Figure 88 Boiling Test Sample M9, 150-Mesh Sintered Nickel Screen 


b. The epoxy bond between the heater and sample maintained good 
thermal contact. 

c * The epoxy edge prevented boiling heat transfer losses from the 
edges of the sample backup plate and heater. Edge boiling could 
also inhibit the flow of liquid to the wick surface. The edge of 
the wick was not covered with epoxy. 

d. The Transite prevented boiling from the bottom of the heater. The 
Transite was held in place by the epoxy on the edges of the assembly 
rather than by epoxy on the underside of the heater. The epoxy 
could have decomposed if used on this surface, since the bottom 
surface of the heater was at the highest temperature of any surface 
in the boiling apparatus. 

3. Instrumentation 

Six chromel-alumel thermocouples were located in the underside of the 
sample backup plate. In addition two thermocouples were located in the 
vapor region above the sample, and two immersed in the liquid. Also, 
two sets of thermocouples were located at different planes in the Transite 
block on top of the guard heater to enable the control of this heater. For 
each test point, the guard heater was adjusted until the two sets of thermo- 
couples read approximately the same, thus minimizing heat losses from 
the underside of the heater through the Textolite platform. Figure 89 
shows the location of thermocouples installed on each sample and of the 
thermocouples used to regulate the guard heater. Figure 88 shows six 
thermocouples installed on a 150-mesh sintered screen sample. The six 
thermocouple wires installed on the sample were enclosed in sheaths of 
0. 020-inch outside diameter which were welded in 0. 020-inch grooves on 
the underside of the sample backup plate. The thermocouple wires (0. 0025 
inch in diameter) were welded to the bottom of the grooves, and the grooves 
filled with epoxy. Hence, the temperature being measured was at a dis- 
tance of 0. 002 or 0. 005 inch from the liquid-wick- solid interface, since 
the sample backup plates were either 0. 022 or 0. 025 inch thick. 

Other instrumentation included a voltmeter to measure the voltage drop 
across the sample heater, and two ammeters to measure the current through 
the sample and guard heaters. 

Table 14 presents a summary of the instrumentation used in the boiling rig 
with the maximum error value of each instrument. 

4. Power Supplies 

Power for the boiling rig is supplied by a 12 KVA rectifier with a three- 
phase 480-volt input and a controllable DC output of up to 750 amperes at 
12 volts, or 1500 amperes at 6 volts. Power for both the guard heater 
and the heater located between the primary and auxiliary tanks is controlled 
by 120- volt input powerstats. The secondary tank heater is a 5 KW 230- 
volt resistance type, controlled by a powerstat. 
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dimensions in inches 


Figure 89 Location of Thermocouples for Sample and Guard Heater 
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TABLE 14 

Boiling Test Instrumentation 


Instrument 

Number and 
Measured Quantity 

Location 

Range 

of 

Meter 

Range 

of 

Tests 

Maximum 

Error 

voltmeter 

1, heater voltage, E 

across heater 

0-7. 5v 

0. 5-4v 

0. 02 5v 

ammeter 

1, heater current, I 

before heater 

0-1000 amps 

75-600 amps 

5 amps 

ammeter 

1, guard heater 
current , Ig 

before guard 
heater 

0-30 amps 

2-10 amps 

0. 25 amps 

C-A thermocouple 

4, guard heater, T-^, 

Trk 

(see Figure 92) 

0 - 800°F 

100-650°F 

1. 5°F 


block temp. T 3 , T^ 



- 


C-A thermocouple 

6, sample temp. T 5 , 

t 6 , t 7 , t 8 , t 9 , t 10 

(see Figure 92) 

0 - 800°F 

1 00-65 0°F 

1. 5°F 

C-A thermocouple 

2, primary inlet 
temp. T^, T ^2 

inlet tube to 
primary tank 
approx. 1 inch 
from tank wall 

0 - 800°F 

100-200°F 

1. 5°F 

C-A thermocouple 

2, vapor temp. T ls , 
T 14 

1/2 inch above 0 - 800°F 
sample, 1/2 inch 
from either end, 
half way across in 
width direction 

200-220°F 

1. 5°F 

C-A thermocouple 

2, primary tank 
water temp. T^g, T^g 

adjacent to 0 - 800°F 

Textolite plat- 
forms, 1/2 inch 
from sample ends, 
immersed in water 
3/4 inch 

100-200°F 

1. 5°F 

C-A thermocouple 

1, busbar temp. 
T 

17 

probe for busbar 
temperature at 
any location 

0 - 800°F 

100-650°F 

1. 5°F 


D. Test Procedure 


The boiling tests were started by filling the secondary tank with the test fluid to 
the approximate level desired in the primary tank. The hose clamp in the line 
connecting the two tanks was opened, allowing the fluid to flow from the secondary 
tank through the heater to the primary tank. The heater was adjusted to fill the 
primary tank with fluid at the desired temperature. Fine adjustments of the 
liquid level in the primary tank to change the degree of liquid saturation in the 
wick sample could be made by altering the adjustable feet on the primary tank 
or the level of liquid in the secondary tank. 

Power was supplied to the primary heater, and the input varied to obtain the 
desired heat fluxes at which data was to be taken. At each test point the guard 
heater was adjusted by means of a 120V powerstat until the temperatures of the 
top two Transite-block thermocouples read approximately the same as the bottom 
two (indicating neglible heat loss from the bottom of the heater through the 
Textolite block). The temperatures were recorded after stabilizing to a point 
where there was no change over a three-minute period. The heater voltage 
and current and the guard heater current were also recorded. Adjustments 
in the liquid level could be made between data points, if necessary. (The 
decrease in water level due to boil-off was 0. 018 inch for a 10-minute run at 
a heat flux of 300, 000 BtuAr ft 2 . ) 

As power input was increased and sample temperature rose above the saturation 
temperature, vapor bubbles formed at nucleation sites within the wick. Further 
increase in sample temperature resulted in the formation of larger and more 
numerous bubbles until for some cases the sample was finally blanketed with 
vapor, indicating film boiling conditions. At the onset of film boiling, as evi- 
denced by a sharp rise in sample temperature, the test was terminated. In 
other cases, the test was terminated when the temperature limitation (about 
450°F) of the epoxy bond was indicated by any guard heater thermocouple. 

E. Calculation Procedure for Heat Flux and AT sa t 


The heat flux of the heater was computed from the measured current and voltage 
drop across the sample, and the area of the sample. In order to determine the 
sample heat flux, the heat flux value of the heater was adjusted for heat loss by 
using a calibration curve for rig heat loss. 

The wall temperature at the liquid-wick-solid interface can be calculated using 
the equation 


Twali 


T - 

m 


Q Ax 
A *8 


( 102 ) 
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where T 

<A 


Ax 


*8 


measured temperature of backup plate 
sample heat flux 

backup plate thickness from thermocouple locations to the 
liquid- wick-solid interface (0. 002 or 0. 005 inch depending 
upon the sample) 

thermal conductivity of the backup plate 


The temperature at the liquid-wick-solid interface minus the fluid saturation 
temperature, AT sa ^, is calculated as follows: 


AT sat = Twal l - T sat < 103 > 

where T sa ^. is the fluid saturation temperature. 

The test data of AT g . vs q presented in the next section is based on averages 
of T wa ij calculated from the several measured temperatures. 

F. Test Results 


The following sections contain a discussion of the test results obtained from the 
study on boiling from wick covered surfaces. The effects of liquid level (degree 
of submergence), orientation, surface modification by grooving the wick sur- 
face, fluid properties and wick structure are shown by comparing mean lines 
drawn through the data points for each set of data. All tests were run at atmos- 
pheric pressure. Boiling heat transfer data was obtained for an Inconel flat 
plate in the horizontal position in order to have a reference with which to com- 
pare the results of wick boiling tests. Flat plate data is shown in Figure 90, a 
plot of heat flux versus ATsat- Table 15 presents an outlined description of the 
boiling tests. 

1 . Effect of Degree of Submergence for Wick in Horizontal Position 

The effect of water level, or degree of submergence, is demonstrated in 
Figure 91 for a sintered 150-mesh nickel screen tested in a horizontal 
position. The degree of submergence is the ratio of the thickness of wick 
that is submerged to the total thickness of the wick. The boiling character- 
istics of the sample did not change significantly as the degree of submergence 
changed from 100 per cent to 25 per cent. As shown in Figure 91, the two 
curves representing 100 per cent submergence are not similar. The dif- 
ference is due to the fact that this sample was rebonded to the heater between 
Tests No. 1 and 5, and indicates that boiling heat transfer characteristics 
of a wick covered surface depend strongly on the physical handling of the 
sample. During these tests, it was determined, by visual observation, 
that all vapor bubbles emanated from the edge surface with none going 
through the top surface. In the case of 100 per cent submergence the liquid 
could be seen on the top and no bubbles were observed emanating from the top. 
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4 T^at = T W all“T S at( 0 F) 

Figure 90 Heat Flux vs AT sat for Horizontal Flat Plate 
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TABLE 15 


Summary of Boiling Tests 






Orientation 


Range of 

Maximum Thermocouple Temperature 


Results Shown 


Test 

(Angle to 

Degree of 

Fluid Bulk 

Variation From Mean Reading 


Test No. In Figure No. 
0 94, 110 

Sample 

Fluid 

(Horizontal - Deg. ) 

Submergence (%) 

Temo 

AT at a. AT - °F and o - Btu/Hr Ft 2 

Comments 

flat plate 

water 

0 

100 

q <262.000 

2.5 at 224.000 

uniform boiling on surface 





T * 140*193 



1 95 

sintered 150-mesh 
nickel screen - M9 



0 

100 

q < 89.000 

23 at 15, 000; 42 at 69, 000; 27 at 89, 000 

all vapor evolution through sides, non-uniform 


T * 145*196 temperature distribution due to preferred paths 

of vapor escape 


2 


9G 


24 


5/8” of wick q <48,000 2 at 1 D, 000: 5 at 48, 000 for TC's at sub- 
end submerged T » 111-*-! 51 merged end - only one TC used at other end 


3 96, 111 

4 96 


45 


90 


1/4” of wick q <41,000 

end submerged T ~ 97*144 


3/8” of wick q < 33.000 

end submerged T ~ 98*168 


liquid appeared to be oscillating (period 1 sec) 

0 at 16,000; 5 at 41,000 for TC's at sub- liquid near the wick surface one instant and 

merged end - only one TC used at other end away from it the next - no vapor bubbles were 

observed 

1 at 10,000; 5 at 33,000 for TC's at sub- 
merged end - only one TC used at other end 


95, 110 


0 


100 


q < 43,000 10 at 11, 000; 45 at 43,000 

T « 126*166 


6 


95 


0 


50 


q <35,000 14 at 9,000; 47 at 28,000 

T « 112*132 


same as test no. 1 


7 


95 


0 


25 


q <43,000 3 at 10,000; 47 at 3G.OOO; 34 at 43,000 

T ~ 134*145 


8 97,110 sintered 50-mesh 

nickel screen - M7 


0 


100 


q <90,000 5 at 16,000; 12 at 90,000 

T ~ 151*189 


same as test no. 1 


9 

97 

same as test no. 8 
with 5 layers of screen 
removed 

0 

100 

q < 58, 000 
T ~ 96*180 

2 at 11,000; 25 at 58, 000 vapor evolution through top surface at discrete 

sites - some through sides 

10 

98, 100 
109, 110 

sintered nickel 
powder - M2 

0 

100 

q < 129,000 
T ~ 118*185 

— 4 at < 80, 000; 18 at 117, 000; 43 at 129, 000 vapor evolution predominantly through top sur- 
face, some through sides - due to nonuniformities 
in wick structure. One end of sample had 

11 

98 

1 

0 

50 

q < 93,000 
T« 121*181 

2 at 22,000; 5 at 67,000; 29 at 93,000 higher limiting heat flux than other end, as shown 

In figure 9 3for tests no. 11 & 12. Not shown for 
test No. 10 since does not change trend of curve 


appreciably 
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TABLE 15 (Cont'd) 


Test No - 

12 

13 

14 

15 
1G 

17 

18 

19 

20 
21 
22 

23 

24 


Results Shown 


Test 

In Figure No. 

Sample 

Fluid 

98 

sintered nickel 

water 


powder - M2 

1 


98 

99,111 

99 

100 

101 , 108 , 

110 

101 

101 
103 

103 

104, 110 sintered stainless steel 

fiber - HU 

104 I 


sintered nickel 
powder - M2 - with 2 
grooves lengthwise 
along face 

sintered nickel 
fiber - H6 

sintered nickel fiber 
HG with 2 grooves cut 
lengthwise along face 


Orientation 
(Angle to) 

(Horizontal - Deg. ) 

Degree of 
Submergence (%) 

Range of 
Fluid Bulk 
Temp 

Maximum Thermocouple Temperature 
Variation From Mean Reading 
AT at q, AT - °F and q - Btu/Hr-Ft 2 

Comments 

0 

25 

q < 93, 000 
T w 133 -*-178 

2 at 22,000; 8 at 79,000; 56 at 93,000 

same as tests no. 10 & 11 

90 

100 

q < 77,000 
T » 131-^-212 
q >77,000 
T = 212 

1 at 23, 000; 8 at 77,000; 13 at 93,000 

same as test no. 10 through 12 

45 

1/4" of wick 

q < 74,000 

1 at 41,000; 4 at 74,000 - TC's at submerged 

fluid motion in the wick that decreased as the 


end submerged 

T« 99 -*-131 

end; 1 at 41, 000; 5 at 74, 000 - TC's at other 
end 

distance from the tank liquid level increased 
was observed - could not determine whether 
motion was due to an oscillating interface or 

90 

l/2" of wick 
end submerged 

q < 39,000 
T~ 94~»-133 

0 at 19,000; 5 at 39, 000 - TC's at submerged 
end; only 1 TC used at other end 

vapor bubbles 

0 

100 

q < 180,000 
T ~ 95 ••■191 

3 at 30,000; 11 at 180,000 

q of 10, 000 - 20, 000 vapor evolution from wick 
surface and grooves - as q increased vapor 
evolution predominantly from grooves - at max 
q (180, 000) sample nearly obscured due to 
vigorous boiling 

0 

100 

q < 103, 000 
T ~ 123-^193 

1 at 21,000; 3 at 103,000 

vapor evolution through top surface ■ 

0 

100 

q < 180, 000 
T =*> 104 -*-182 

1 at 15,000; 6 at 180,000 

vapor evolution predominantly from grooves - 
some from sides 

0 

50 

q < 128,000 
T » 1 19-^184 

2 at 30,000; 7 at 128,000 

vapor evolution from grooves and sides 

45 

1/2” of wick 

q < 20, 000 

2 at 7, 000; 2 at 20, 000 - TC's at submerged 



end submerged 

T 3; 99 ♦117 

end; only 1 TC used at other end 

same as tests no. 14 & 15 

90 

1/2" of wick 

q < 21,000 

1 at 14,000; 4 at 21,000 - TC's at submerged 



end submerged 

T ~ 114 *-135 

end; 1 at 14,000; 0 at 21,000 TC's at other end 

0 

100 

q < 141.000 
T « 111 -*-185 

1 or 2 at 105, 000; 10 at 141, 000 


0 

50 

q < 58, 000 
T» 103 -►171 

1 at 17, 000; G at 58, 000 

vapor evolution predominantly through top 
surface - some through sides 

0 

25 

q < 71,000 
T « 98-M74 

<2 at <71,000 
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Test No. 
25 

2G 

27 

28 

29 

30 

31 

32 

33 

34 

35 

36 

37 


TABLE 15 (Cont'd) 







Orientation 


Range of 

Maximum Thermocouple Temperature 


Results Shown 



Test 

(Angle to) 

Degree of 

Fluid Bulk 

Variation From Mean Reading 


In Figure No. 

Sample 

Fluid 

(Horizontal - Deg. ) 

Submergence (%) 

Temp 

AT at q, AT - "F and q - Btu/Hr-Ft 2 

Comments 

105 

sintered stainless steel 

water 

0 

100 

q < 74,000 

1 at 11,000; 2 at 38,000; 2 at 74,000 

at q = 49,000 all vapor evolution was througn 


fiber - 

H13 





T« 118—183 


top surface; at q = 74,000 some vapor also 
evolved through sides 

105 





0 

50 

q < 49, 000 

2 at 11,000; 3 at 49,000 

at q = 20,000 all vapor evolution was through 








T =s 115—171 


top surface; at q = 28, 000 some vapor also 
evolved through sides 

105 





0 

25 

q < 78,000 

1 at 11,000; 4 at 78,000 

vapor evolution through top surface - 








T « 120—177 


possibly some through sides 

105, 106 





0 

100 

q < 132,000 

2 at 23,000; 19 at 132,000 

vapor evolution predominantly through top 

109, 110 







T» 140—195 


surface - some at sides 

106 





45 

100 

q < 70,000 
T = 104—210 

3 at. 70,000; 6 at 70,000 

vapor evolution through top surface 

106 





90 

100 

q < 32,000 
T» 95—149 

1 at 12,000; 4 at 32,000 

same as test no. 28 

107, 111 





45 

1" of wick end 

q < 74,000 

only 1 TC used at submerged end; 4 at 








submerged 

T as 100—141 

18,000; 1 at 74,000 - TC's at other end 

same as test nos. 14 and 15 

107 





90 

1” of wick end 

q < 27, 000 

only 1 TC used at submerged end; 2 at 








submerged 

T as 97—139 

17,000; 3 at 27,000 - TC's at other end 


108 

smooth Inconel flat 

Freon 113 0 

100 

q < 44,000 

4 at 18,000; 3 at 114,000 

boiling over entire surface 


plate 






T a: 78—117 
q > 44,000 
T = 117 



108 

roughtened Inconel 
flat plate 



0 

100 

q < 37,000 

4 at 23,000; 4 at 112,000 

boiling over entire surface 








T w 79—117 
q > 37,000 
T = 117 



108, 109 

sintered nickel 



0 

100 

q < 37, 000 

9 at 23, 000; 28 at 37, 000; 120 at 61, 000 

vapor evolution across top surface 


powder 

- M2 





T ss 75—117 
q > 37,000 
T = 117 



108,109 

sintered nickel fiber 



0 

100 

q < 76,000 

2 at 32,000; 8 at 73,000; 22 at 76, 000 

same as test no. 35 


- HG 






T » 81—111 



108,109 

sintered stainless steel 


1 

0 

100 

q < 78. 000 

6 at 19,000; 35 at 51,000; 113 at 78,000 

same as test no. 35 


fiber - H13 76 ♦■111 


id aH/rua- c _oixb 



^"^SAT = TwALL~ T sat (°F) 


Figure 91 Heat Flux vs AT sa t for Sintered 150-Mesh Nickel Screen Wick 
(Sample M9) in Horizontal Position 
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Figures 94, 97, 100 and 101 show curves of heat flux vs AT g for several 
other wick sainples, with degree of submergence as a parameter. These 
curves indicate that for both a sintered metal powder wick (Figure 94) and 
sintered metal fiber wicks (Figures 97, 100, and 101) a decreasing water 
level results in a reduced heat transfer rate for a given AT gat . Hence, 
the boiling heat transfer coefficient for these samples decreases as the 
degree of submergence decreases. During these tests, nucleation sites 
were observed scattered across the surface of the wick and vapor evolution 
occurred at the face of the sample, rather than at the edges, as occurred 
with the sintered 150-mesh nickel screen. 

Some of the data shown in Figure 97 was obtained using a sample with 2 
grooves cut lengthwise along the face. A sketch of this sample is shown 
in Figure 98. The presence of the grooves did not alter the result that 
heat flux at a particular AT gat decreased as degree of submergence de- 
creased. 

2. Effect of Thickness of a Sintered Metal Screen Wick 


The effect of wick thickness for a sintered metal screen was determined in 
two tests on a sintered 50-megh nickel screen sample. These were a 
horizontal, 100 per cent submergence test on a 0. 106-inch thick sample 
made of 7 screen layers, and a horizontal 100 per cent submergence test 
on the same sample with the top 5 screen layers removed. During the 
first test, no bubbles were observed leaving the top surface of the sample. 
Instead, as in the 150-mesh screen sample test, the vapor bubbles left 
the wick only at preferred sites along the edges of the sample. On the 
other hand, in the test with only 2 layers of screen, vapor bubbles were 
seen escaping from the top surface of the sample at discrete sites. The 
results of these tests are shown in Figure 93. The boiling heat transfer 
coefficient of the 2-layer thick sample is significantly higher than that of 
the 7-layer thick sample, probably due to the ease of vapor escape rather 
than to any difference in nucleation within the wick. 

3. Effect of Inclination With End of Wick Submerged 

The object of these tests was to determine if the orientation of gravitational 
forces affect the boiling characteristics of a wick covered surface. To 
accomplish this , tests were made in which a sample was inclined to the 
horizontal, with a section of the wick, from 1/4 to 1 inch, submerged in 
the liquid. The liquid wicked up the sample due to capillary forces and 
was vaporized. The results of these tests are shown in Figures 92, 95 , 

99 and 103, for four test samples. 
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Figure 92 Heat Flux vs AT sa t for Sintered 150-Mesh Nickel Screen Wick 
(Sample M9) at Angles to Horizontal 
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Figure 93 Heat Flux vs AT sat for Sintered 50-Mesh Nickel Screen Wick 
(Sample M7) 
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o 

10 
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II 
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SHADED SYMBOL - POINT OF DECREASING HEAT FLUX 
] 


AVERAGE OF OTHER 
THERMOCOUPLES FOR 
SAME TEST, 


TEST FLUID -WATER 
PRESSURE -ATMOSPHERIC 



AVERAGE OF 2 THERMO- 
COUPLES LOCATED AT 
ONE END OF SAMPLE 
(SEE COMMENTS 
IN TABLE 15) 


O 40 80 

^ Tsat = T wall - T sat (°F) 
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Figure 94 Heat Flux vs AT sat for Sintered Nickel Powder Wick (Sample M2) 
in Horizontal Position and Completely Submerged in Vertical 
Position 
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Figure 95 Heat Flux vs AT sa t for Sintered Nickel Powder Wick (Sample M 2 ) 
with Partial Submergence at Angles to Horizontal 
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Figure 96 Heat Flux vs AT sa t for Grooved and Ungrooved Sintered Nickel 
Powder Wicks (Sample M2) 
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sat = T W all- "I"sat (°F) 


Figure 97 Heat Flux vs AT sa f. for Sintered Nickel Fiber Wick (Sample H6) 
in Horizontal Position 
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section 


unheated end sections are 
crushed under clamping 
blocks (see Figure 85) and 
do not function as wicks 


groove depth 1/16 inch 

nominal sample thickness 0. 10 inch 


all dimensions in inches 


Figure 98 Dimensions of Grooved Samples 



“100 “80 -60 -40 -20 0 20 

4Tsat = Twall~ Tsat (°F) 

Figure 99 Heat Flux vs AT sa t for Grooved Sintered Nickel Fiber Wick 
(Sample H6) at Angles to Horizontal 
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Figure 100 Heat Flux vs AT gat for Sintered Stainless Steel Fiber Wick 
(Sample Hll) 
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Figure 101 Heat Flux vs AT sa £ for Sintered Stainless Steel Fiber Wick 
(Sample H13) in Horizontal Position 
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^T sat = Twall “ "T S at (°F) 

Figure 102 Heat Flux vs AT sa t for Sintered Stainless Steel Fiber Wick 
(Sample H13) Completely Submerged at Angles to Horizontal 
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^Tsat =T wall -T sat (°F) 


Figure 103 Heat Flux vs AT sa t f° r Sintered Stainless Steel Fiber Wick 
(Sample H13) Partially Submerged at Angles to Horizontal 


PAGE NO. 161 




q XlO BTU/HR FT 


280 


<M 


240 


200 


i 



TEST 

NO. 

0 

33 

V 

34 

A 

35 

o 

36 

O 

37 


SAMPLE | | 

- SMOOTH FLAT PLATE 

- ROUGHENED FLAT PLATE j 
M2 - SINTERED NICKEL POWDER 

H 6 - SINTERED NICKEL FIBER 
H 13 -SINTERED STAINLESS STEEL FIBER 


DEGREE OF SUBMERGENCE - 100% 
ORIENTATION - HORIZONTAL 


UNSHADED SYMBOL -POINT OF INCREASING HEAT FLUX 
SHADED SYMBOL -POINT OF DECREASING HEAT FLUX 



Figure 104 Results of Wick Boiling Tests with Test Fluid Freon 113 


PAGE NO. 162 



Each figure shows a plot of heat flux, q, versus AT ga |. based on sample 
temperatures calculated at two locations; from an average of the thermo- 
couples located closest to the surface of the test fluid (Thermocouples No. 

5 and 8 - Figure 89), and an average of the thermocouples located farthest 
from the surface of the test fluid (Thermocouples No. 7 and 10 - Figure 
89). In some cases, the lower thermocouples were located below the liquid 
surface. 

Figure 92 is a graph of heat flux vs AT sa t f° r sintered 150-mesh nickel 
screen. It is shown that AT sa ^ at the location farthest from the liquid 
surface is greater than that at the location nearest the liquid surface for a 
given average heat flux. The temperature variation is probably due to 
gravity forces and escaping vapor opposing the flow of fluid in the wick. 

The test fluid has to flow by capillary wieking about an additional 2 1/4 inches 
beyond the lower thermocouples to reach the thermocouples located farthest 
from the test fluid. It is likely that only the smallest pores near the top 
of the wick are filled with liquid, resulting in the wick containing considerable 
vapor at this point. This causes a higher resistance to heat transfer, and 
thus a higher sample temperature. A possible contributing cause for capil- 
lary pumping failure of the large pores in the upper region of the wick is 
that vapor generated in the lower regions probably has difficulty in escaping 
from the wick, thus causing a greater resistance for liquid flow to the upper 
region. This is somewhat substantiated by the observation in the horizontal- 
orientation fully-submerged test on the same sample (Test No. 1) that the 
vapor only escaped through the wick sides. Gravity effects may also be seen 
by comparing the temperature of the sample at the end farthest from the 
water for the three tests shown in Figure 92. This temperature can be 
seen to become progressively higher as the angle is increased and the com- 
ponent of gravity force opposing the flow increased. 

Figures 95, 99, and 103 show curves of heat flux versus AT sat for three 
different samples inclined both 45 degrees and 90 degrees to the horizontal 
with one end submerged in the water tank. These curves are for sintered 
nickel powder, sintered nickel fiber and sintered stainless steel fiber wicks, 
respectively. In all cases the average thermocouple reading at the top of 
the sample was higher than at the lower end of the sample. However, no 
appreciable difference was noticed between the 45-degree angle test and 
the 90-degree angle test, except for the sintered nickel powder sample at 
heat fluxes greater than 25, 000 Btu/hr ft 2 (see Figure 95). It should be noted 
that for the other two samples there is no data at the two different angles for 
heat fluxes greater than 25, 000 Btu/hr ft 2 . The results of these tests can be 
explained in the same manner as those of the sintered 150-mesh nickel screen. 
That is, the additional capillary force needed to overcome gravity effects causes 


PAGE NO. 163 


some of the larger pores to fail as capillary pumps. The lower limiting 
heat flux of sample M2 in the 90-degree inclined case as compared to the 
45-degree inclined case was probably due to increased effects of gravity 
at higher heat fluxes , because of the greater requirements of the wick to 
overcome the frictional pressure loss. The frictional pressure loss effects 
are compounded at a higher heat flux since both greater amounts of vapor 
exist and a higher liquid flow is required. 

For all of the samples tested at an inclination, the tests on the same sam- 
ples in the fully submerged horizontal position indicated lower AT g ^ at 
the higher heat fluxes than the average of the upper thermocouples^ d for 
the inclined cases at the same heat flux. This can be expected, since for 
the horizontal cases the liquid has to wick, at most, one-half of the sample 
width (5/16 inch), rather than a large fraction of sample length (2 5/8 inches). 

4. Effect of Inclination With Wick Fully Submerged 

Figure 94 contains the results of a test on a sintered nickel powder wick 
(M2) in a vertical position and completely submerged. This curve can be 
compared to the curve contained in the same figure with the sample in a 
horizontal position and completely submerged. The latter case produced 
both a lower AT sa t f° r a given heat flux in the boiling regime and a higher 
limiting heat flux. Figure 102 shows a similar effect of inclination on 
AT sa t a t a given heat flux for a sintered stainless steel fiber wick when 
fully submerged. The most likely explanation for these results is that the 
vapor can more easily escape from the wick surface in a horizontal position. 

5. Effect of Grooves Cut in Wick 


Figures 96 and 97 are graphs of heat flux versus AT sa t for both grooved 
and ungrooved wicks oriented horizontally. Figure 96 shows the effect of 
grooving a sintered nickel powder wick, and Figure 97 the effect of groov- 
ing a sintered nickel fiber wick. In both cases it can be seen that grooving 
the wick both reduces the degree of wall superheat required for a given 
heat flux and increases the limiting heat flux. As shown in Figure 96 the 
limiting heat flux for an ungrooved sintered nickel powder wick is about 
130, 000 Btu/hr ft 2 , while that for a grooved sample is greater than 
180, 000 Btu/hr ft 2 . The grooved wicks have 2 slots cut lengthwise along 
the face which measure 2 7/8 inch long, 3/64 inch wide and 1/16 inch deep. 
Figure 98 is a drawing of a grooved wick, showing the location of the grooves. 
The improved performance of the grooved samples was probably due to 
smaller resistance for vapor escape offered by the grooves. 
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6. Effect of Fluid Properties 


To determine the effect of fluid properties on boiling from a wick-covered 
surface, tests were made with several samples in a horizontal position at 
100 per cent submergence with Freon 113 as the test fluid. The results of 
these tests are shown in Figure 104, and are replotted in Figure 105 to 
allow for a comparison of tests on the same sample with both water and 
Freon 113 as test fluids. The results of the tests shown on Figure 105 
demonstrate that invariably a higher heat flux was obtained at a given 
AT S at when water was used as the test fluid. The limiting heat flux was 
also significantly higher when water was used as the test fluid as shown 


in the table below. 



Sample 

Test Fluid 
Water 

Test Fluid 
Freon 113 

sintered nickel 
powder wick (M2) 

128, 000 Btu/hr ft 2 

63,000 

sintered nickel 
fiber wick (H6) 

>104,000 

75, 000 

sintered stainless 
steel fiber wick (H13) 

>132, 000 

80,000 


As a means of analytically comparing the values of limiting heat flux 
found experimentally with water and Freon 113 as the test fluids, the fol- 
lowing critical heat flux correlation, proposed by Rohsenow and Griffiths 15 , 
was used: 


^ critical 
p v h VL 


= 143 g 


1/4 




where q - 

P L = 
P v = 



heat flux 

liquid density 

vapor density 

heat of vaporization 

acceleration due to gravity 


(104) 


While this correlation is for a bare surface, rather than a wick -covered 
surface, it allows a means of comparison of critical heat flux for various 
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fluids based on physical properties. Assuming vapor density to be small 
compared to liquid density, from Equation (104) the critical heat flux for 
Ereon 113 divided by the critical heat flux for water is 


w 


critical 
Freon 113 


W cri 
wal 


critical 

water 


VL Freon 113 
l VL water 


(~ 


P \ 6 / \ 

L Freon 113\ / p v Freon 113 \ 


water 


which reduces to 


( \ critical 
q ) Freon 113 

\ critical 
q ) water 


( 




0.4 


p v water 


7 


(105) 


(106) 


when saturation properties at atmospheric pressure are substituted. 


Although the ratio of peak heat fluxes obtained from the data was greater 
than this ratio, the analytical prediction that Freon 113 has a lower criti- 
cal heat flux was verified. 


7. Effect of Wick Structure 


In addition to determining the boiling heat transfer characteristics of a 
wick- covered surface, it is further necessary to evaluate and classify the 
performance of wicking materials. To accomplish this, comparisons were 
made between the boiling heat transfer characteristics of a flat plate, 
sintered screen wicks, a sintered metal powder wick, and sintered metal 
fiber wicks. 

A graph of heat flux versus AT sa £ for several samples tested horizontally 
at 100 per cent submergence with Freon 113 as the test fluid is presented 
in Figure 104. A comparison of the curves shows that the sintered metal 
fiber wicks (H13 and H6) as boiling surfaces are better than the sintered 
metal powder wick (M2), but not as good as smooth or roughened flat plates, 
when the distance the liquid has to travel is insignificant. 

Figure 106 allows the same type comparison as Figure 104, only with water 
used as the test fluid. Included in Figure 106 are results of tests on three 
different sintered fiber wicks, one sintered powder wick, and two different 
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Figure 105 Comparison of Wick Boiling Characteristics of Three Samples 
Tested in Water and Freon 113 
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Figure 106 Results of Wick Boiling Test with Water as Test Fluid 
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Figure 107 Comparison of Sample Boiling Characteristics When Fluid Wicks 
up Sample 


screen wicks. The three sintered metal fiber wicks have the best boiling 
characteristics, followed by the sintered metal powder wick, and then the 
two sintered screen wicks. With water as the test fluid, two of the sintered 
metal fiber wicks (H13 and Hll) perform as well as the flat plate (see 
dashed line Figure 106), contrary to the results of the tests with Freon 113 
as the test fluid and H13 as the test sample. Probably Hll and H13 would 
have shown poorer performance than the flat plate when tested in water if 
a higher heat flux had been obtained. The poorer performance of H13 com- 
pared to the flat plate when tested in Freon 113 may be attributed to vapor 
blanketing of the wick sample (i. e. limiting heat flux being attained) at 
lower heat fluxes than with water. 


In order to compare the boiling characteristics of several wick materials 
when the test fluid has to be wicked through the length of the sample, some 
of the results of tests with the sample inclined at an angle of 45 degrees to 
the horizontal, with a section submerged in the water, were replotted in 
Figure 107. This figure presents curves of heat flux versus AT ga .£, where 
the sample temperatures used in the calculation of AT sa t are averages of 
the temperatures measured by the thermocouples located farthest from the 
liquid surface (Nos. 7 and 10 - Figure 93). Results are shown for the 
sintered 150-mesh nickel screen wick (M9), a sintered nickel powder wick 
(M2), and a sintered stainless steel fiber wick (H13). This figure shows 
that the sintered stainless steel fiber wick had the best boiling characteristics, 
such as a low AT sa t (18°F) for a given heat flux (74, 000 Btu/hr ft 2 ) and a 
relatively high limiting heat flux ( > 74, 000 Btu/hr ft 2 ). The sintered nickel 
powder sample had fairly good boiling heat transfer characteristics (i. e. 
low ATg a t) up to a heat flux of about 62, 000 Btu/hr ft 2 . The limiting heat 
flux for this sample was about 74, 000 Btu/hr ft 2 . The sintered 150-mesh 
nickel screen exhibited the poorest performance, with a very high degree 
of wall superheat (144°F) for a low heat flux (41, 500 Btu/hr ft 2 ). The 
AT sa t for this sample increased gradually as heat flux increased, with no 
sudden break in the curve. This was probably due to the inability of the 
vapor to escape readily, as was previously discussed. 


In order to see if the limiting heat fluxes in these inclined tests were as- 
sociated with failure of the wicks as capillary pumps, Equation (50) was 
rewritten as follows: 



qu x * 
E 


\ 


T « P T S„ 
L o 


K 


P_ sin 0 x ( 10 ?) 

L E 


where, in this case, q' and x E are the evaporator heat flux and length, 
respectively. 
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The term on the left is the capillary pump capability and the terms on the 
right are the frictional loss and gravitational head loss, respectively. 

The experimentally determined values of f> m from the wick rise tests, and 
Kp from the permeability tests, listed below, and the other known values 


were substituted in the above equation to solve for the predicted limiting 

heat flux. 

Sample 

K X (ff 2 ) 

Z m (inch) 

sintered 150-mesh nickel screen wick (M9) 

9. 0 x 10 8 * 

6.9 

sintered-nickel powder wick (M2) 

3.4 x 10 8 

9.7 

sintered stainless steel fiber wick (H13) 

0. 8 x 10 8 

5.3 


* estimated value based on interpolation of values for samples M8 and M10 


The results of the calculation are presented below: 

Predicted Value of 

Sample Limiting Heat Flux-Btu/hr ft 

sintered 150-mesh nickel screen wick 0.2 x 10 6 

g 

sintered-nickel powder wick 1.2 x 10° 

sintered stainless steel fiber wick 4. 6 x 10 8 

The only sample that experimentally exhibited a definite limiting heat flux 
was the sintered nickel powder wick, and this value was considerably lower 
than the calculated value. However, the above analysis assumed liquid 
flowed through all of the pores, whereas in the boiling tests, vapor was 
probably present in a considerable fraction of the pores. If considerable 
vapor were present, the frictional pressure loss would be much higher due 
both to a reduction in flow cross-sectional area and perhaps an increase in 
if the Darcy regime were exceeded. 

It appears from inspecting Figure 107 that the sintered fiber sample would 
attain the highest limiting heat flux, followed by the sintered powder sample. 
The sintered screen sample did not show a distinct limiting heat flux in the 
range tested. However, it appears that it would fail before the sintered powder 
sample. The results of the analysis show the same order in limiting heat 
flux as that mentioned above. This agreement might only be fortuitous since 
it is difficult to ascertain the fraction of pores containing vapor in each of the 
different samples and the effects of boiling on limiting capillary pumping heat 
flux. 


G. Summary and Conclusions 

The test results obtained from a study on boiling from a wick-covered surface 
indicate that: 
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1. Of the porous materials tested, the sintered metal fiber wicks exhibited 
the best boiling characteristics in both a horizontal and inclined position. 
Although the sintered metal screen wicks had the poorest performance, 
they could be used as a wicking surface at relatively low heat fluxes. How- 
ever, even at low heat fluxes a relatively high AT sa t is required. This 
would be undesirable since either a higher system temperature would be 
required for the same heat rejection rate and heat transfer area, or a 
greater heat transfer area would be required for a given heat rejection 
rate and operating temperature. 

2. Vapor in a wick probably causes a significant resistance to the flow of liquid 
through the wick. For a wick with only one end submerged in liquid, the 
greater the distance from the liquid surface, the poorer the boiling per- 
formance. This can probably be attributed partially to vapor blocking liq- 
uid flow up the wick. For fully submerged 50-mesh screen wicks, one 
seven-layer thick screen and the other two layers thick, the thinner wick 
showed a lower AT sa ^ for a given heat flux when both were in horizontal 
position. This was definitely attributed to the inability of the vapor to es- 
cape normal to the surface in the seven-layer case. A further illustration 
of how the ease of vapor escape affects AT sa t was observed in comparing 
tests on two different samples fully submerged, where the samples were 
oriented both horizontally and at angles to the horizontal. In this comparison, 
the horizontal test produced better results which can probably be attributed 

to easier vapor escape. 

3. The greater the degree of submergence in the horizontal position, the 
better the boiling performance for sintered metal powder and sintered 
metal fiber wicks. The boiling test performance of a sintered metal 
screen appeared to be independent of degree of submergence; however, 
in these tests all vapor escaped through the edges of the wick rather 
than through the top surface. 

4. Grooving a wick that has poor boiling characteristics tends to improve 
the performance and increase the limiting heat flux. This is concluded 
from only horizontal testing with full or partial submergence, and further 
testing would be needed to determine whether this conclusion holds for 
other cases. 

5. In the comparison of tests on samples run both with water and Freon 113 

as test fluids, the results showed that much better performance was obtained 
with water. One correlation indicated that critical heat flux is approx- 
imately proportional to vapor density to the 0.4 power, liquid density to 
the 0. 6 power, and the latent heat of vaporization to the 1. 0 power. This 
correlation qualitatively agrees with the test results which indicate that 
water has a higher critical heat flux than Freon 113. 
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IX. CONCLUSIONS 


A. The simplified analysis indicated that the maximum height to which a liquid 
will rise in a vertical wick, j£ m , and the friction factor, Kp are the import- 
ant parameters in determining the capability of a capillary pump in a heat 
pipe. This theory indicates that the larger the value of the capillary pump 
parameter, J?m/Ki, the better the capillary pump, when head effects due to 
gravity are negligible. If an application involves operation at a condition 

in which a capillary pump must work against gravity, consideration must 
be given to both the value of j2m/Kl and ofj^m itself. The important fluid 
properties for capillary pumping were also indicated by the simplified analysis. 

B. The experimentally determined capillary pumping characteristics of the three 
classes of materials indicated that as a class the sintered fiber samples were 
best, followed by the sintered powder samples. The sintered screen samples 
demonstrated the lowest potential as capillary pumps. 

C. Experiments conducted to determine the heat transfer coefficients and limiting 
heat transfer rates for boiling indicate that some of these porous materials 
provide high values of these quantities when adequate paths for vapor escape 
from within the wick are provided. Other samples which do not provide vapor 
removal paths result in values of heat transfer coefficients significantly lower 
than those of a bare flat surface. Three sintered fiber samples were deter- 
mined to have the best heat transfer characteristics of those tested, followed 
by a sintered powder sample and sintered screens. 

D. The analyses and experimental results indicate that vapor-chamber fins 
employing suitable wicks can operate at relatively high evaporator heat 
fluxes and corresponding wick flow rates. 


PAGE NO. 173 



APPENDIX 1 


References 

1. Grover, G. M. , T. P. Cotter, and G. F. Erickson, Structures of Very 
High Thermal Conductance, J. of Appl. Phys. , June 1964, 35, p. 1190 

2. Deverall, J. E. and J. E. Kemme, High Thermal Conductance Devices 
Utilizing the Boiling of Lithium or Silver, Report No. LA-3211, 

April 9, 1965, Los Alamos Scientific Laboratory, Los Alamos, 

New Mexico 

3. Deverall, J. E. and J. E. Kemme, Satellite Heat Pipe, Report No. 

LA-3278- MS, April 20, 1965, Los Alamos Scientific Laboratory, 

Los Alamos, New Mexico 

4. Haller, H. C. , S. Lieblein, and B. G. Lindow, Analysis and Evaluation 
of a Vapor-Chamber Fin-Tube Radiator for High- Power Rankine Cycles, 
NASA TN D-2836, June 1964, NASA Lewis Research Center, Cleveland, 

Ohio 

5. Adamson, Arthur W. , Physical Chemistry of Surfaces, Interscience, 

1960, pp 4-6 

6. Cotter, T. P. , Theory of Heat Pipes, Report No. LA-3246-MS, Los 
Alamos Scientific Laboratory, Los Alamos, New Mexico 

Scheidegger, Adrian E. , The Physics of Flow Through Porous Media, 
Macmillan, 1960, pp 68-90 

8. Fox, H. W. , E. F. Hare, and W. A. Zisman, Wetting Properties of 
Organic Liquids on High Energy Surfaces, J. Phys. Chem, 1955, 59, 
p. 1098 

9. Gorring, R. L. and S. W. Churchill, Thermal Conductivity of Heterogeneous 
Materials, Chemical Engineering Progress, Vol. 57, No. 7, pp. 53-59, 

July 1961 

Carman, P. C. , Flow of Gases Through Porous Media, Butterworths, 

1956, p. 8 

11. Ginwala, K. , T. A. Blatt, and R. W. Bilger, Engineering Study of Vapor- 
Cycle Cooling Components for Space Vehicles, Tech. Doc. Repts. No. 
ASD-TDR-63-582, Sept. 1963, AF Flight Dynamics Laboratory, Air Force 
Systems Command, Wright- Patterson Air Force Base, Ohio, pp 120-148 


PAGE NO. 174 



Allingham, W. D. , and J. A. McEntire, Determination of Boiling Film 
Coefficient for a Heated Horizontal Tube in Water- Saturated Wick 
Material, J. Heat Transfer, Feb. 1961, pp 71-76 

Costello, C. P. , and E. R. Redeker, Boiling Heat Transfer and Maximum 
Heat Flux for a Surface with Coolant Supplied by Capillary Wicking, Chem. 
Engr. Progress Symp. Series, 1963, No. 41, 59, pp 104-113 

Costello, C. P. , and W. J. Frea, The Role of Capillary Wicking and 
Surface Deposits in the Attainment of High Pool Boiling Burnout Heat 
Fluxes, AIChE Journal, May 1964, 10, pp 393-368 

Rohsenow, W. M. , and P. Griffith, Correlation of Maximum Heat Flux 
Data for Boiling of Saturated Liquids, AIChE - ASME Heat Transfer Sympo- 
sium, Louisville, Kentucky, March 1955 


PAGE NO 


175 



APPENDIX 2 


Ap 

Arp 

b 


C L 

D 

D C 

Pk 

D 

F 


AR 


go 

h 

hVL 

J 

k 


Ki 

kb 

^ m 

M 

N 

P 

Pa 

Q 

Qmax 

q 

<*' , 
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Rc 

r 

f 

T 

TM 

T sat 

Twall 

u 

U L 

V 

v T 

w 

w T 

X E 

X M 

x T 


Nomenclature 
Latin- Letter Symbols 

cross-sectional area of pores in wick sample normal to liquid flow 
total cross-sectional area of wick sample normal to liquid flow 
dimension of vapor-chamber fin element in z - direction 
specific heat of an incompressible liquid 
diameter 

equivalent capillary pore diameter defined by Equation (99) 
equivalent frictional pore diameter defined by Equation (98) 
mean pore diameter 
force 

free flow area ratio, defined as Ap/A T 
local acceleration of gravity 
proportionality constant in Newton r s second law 
enthalpy 

latent heat of vaporization 

proportionality constant relating work energy to heat energy 
thermal conductivity 

wick friction factor, defined by Equation (25) 
thermal conductivity of wick backup plate 

maximum height to which a liquid will rise in a vertical wick sample 
momentum PL h VI/ r 

dimensional liquid parameter, — ^ — 
pressure 

atmospheric pressure 
heat energy per unit time 
maximum fin heat rejection rate 
time rate of heat transfer per unit area (heat flux) 
heat flux in evaporator section 
radius of curvature of a liquid-vapor interface 
minimum radius of curvature at the liquid-vapor interface that a 
given wicking material will support 
radius of curvature of the liquid-vapor interface in pore C 
pore radius 
mean pore radius 
temperature 

measured temperature of wick backup plate 

saturation temperature of fluid 

wall temperature at the liquid-wick-solid interface 

velocity 

average liquid flow velocity in the pores of wick 

volume of pores 

total volume of wick 

mass flow rate 

total weight of wick 

total length of evaporator section of vapor-chamber fin 
total length of condensing section of vapor-chamber fin 
total length of condensing and evaporating sections of vapor 
chamber fin 
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Greek- Letter Symbols 


4 


a(r) 

pore size distribution function 

AP 

pressure drop 

^ T sat 

temperature difference between saturated vapor and the liquid- 
wick-solid interface 

Ax 

wick sample length in flow direction 

8 

liquid-wick thickness 

€ 

porosity, defined as total void volume divided by total sample volume 

e 

angle fin makes with the horizontal, positive when heated section 
is at elevated end and negative when heated section is at lower 
end 

<t> 

liquid- solid contact angle measured within liquid 

v 

absolute viscosity 

p 

density 

a 

liquid- vapor surface tension 

Superscripts 

c 

due to convection 

k 

due to conduction 

Subscripts 

B 

backup plate 

d 

friction and inertia 

e 

equivalent 

F 

pertaining to Freon 113 data 

f 

due to friction 

g 

pertaining to gravity 

L 

liquid 

med 

median 

P 

pore, or due to pressure 

S 

solid 

T 

total 

V 

vapor 

w 

pertaining to water data 

WK 

pertaining to conditions of wicking rise tests 

X 

at a given x 

x + dx 

at a given x + dx 

a 

pertaining to capillary effects 

1 

at point 1 

2 

at point 2 
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APPENDIX 3 

Derivation of Fin Equation Including Effects 
of Gravity and Evaporator Section 

The model of the vapor-chamber fin shown in Figure 6 is modified by: 

1) Including gravity forces ( The fin is inclined at an angle 0 measured 
from the horizontal. The angle 0 is defined as positive when the heated 
section is at the elevated end and negative when the heated section is at 
the lower end. ) 

2) Including the evaporator pressure drop by assuming the same theoretical 
model in the evaporator as in the condenser section. 

The model used in the second modification can be considered suitable at a low 
heat flux where heat is conducted through the wick-liquid composite to the 
liquid-vapor interface where evaporation takes place. At a high heat flux the 
boiling process in the wick structure is more complicated, and this simplified 
model may not apply. 

With gravity forces included, a gravity force, Fg, is added to the forces acting 
on the differential element shown in Figure 7. The gravity force is written as 

Fg = - Pl « b 5 sin 0 dx (3-1) 

g o 

and acts in the direction of flow for values of 0 between 180 and 360 degrees, 
and opposite to the direction of flow for values between zero and 180 degrees. 

The differential momentum equation (28) becomes 


- eK^L u L dx _ p L s in 0 dx - -j^-~ (UL) dx (3-2) 


dR 


R 


g 


o 


g Q dx 


The combined energy- continuity equation must be modified to account for the 
rate of work done by the gravity force on a volume element, and Equation (33) 
becomes 

dh du 

u p T — - - (h v - h T ) p — + -p* — 

L L dx E L dx 5c 


P L U L 


sin 0 = 0 


(3-3) 


where J is the proportionality constant relating work energy to heat energy. 
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In the condenser section Equation (39) becomes 
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/g sin 0 


U L = 


{ 6 P. 


Vl ' ' ' l 

If gravity forces are neglected 
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J h, 
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u„ 


dx 


(3-4) 


VL 


u. 


h 8 e p r 
VL L 


x 


or du = 


\ L a 


cp. 


dx 


(3-5) 


from Equation (3-4). For values of 6 between zero and 180 degrees, the 
velocity of the liquid in a gravity field will be less than the value calculated 
from Equation (3-5). The maximum value of the second term on the right side 

of Equation (3-4) divided by the first term is approximately 


x 


which is much less than one. 


g 0 2J h VL 


Therefore, the gravity term in Equation (3-4) may be neglected, and the expres- 
sion for liquid velocity is unchanged when gravity forces are present. 


Equations (3-2) and (3-5) can be combined to give an energy-continuity-momentum 
integral equation, valid in the condenser section 



R o 

x = o 



p sin 0 


dx 



x dx 


(3-6) 


where x is the distance from the origin shown in Figure 6 to any point in the 
condenser section. 


The velocity of the liquid in the fin may be plotted vs fin length as shown in the 
sketch below. 
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Fin Length 

Integration of Equation (3-5) beween x = 0 and x = x m shows that 


q x_ 


u 


m 


max 


h VL ,,P L 


(3-7) 


and integration using the boundary condition 


q x 
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x = x 


m 


h 5 ep 
VL L 


and u = u_ yields 

x - x L 


Q x 


m 


” L = \l" P L (X ' Xm)+ h VL , ' p ] 


(3-8) 


where x is a value greater than x m and q' is the evaporator heat flux. 

Since the heat transferred in the evaporator is equal to the heat transferred in 
the condenser, 


q x 


m 


or q ' 


-q' ( X T “ x m) 
_ q x m 
X T " x m 


(3-9) 
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Substitution of Equation (3-9) into Equation (3-8) gives an expression for the 
liquid velocity in the evaporator section 


q x. 


u. 


m 


Vl 4 • P L 



(3-10) 


Equations (3-2) and (3-10) can be combined to give an an energy continuity-momentum 
integral equation valid in the evaporator section. 
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sin 
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dx 
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m 


d (x T - x) 2 (3-11) 


To evaluate the maximum output of the capillary pump, the liquid radius of 
curvature must be very large at x = 0 (i. e. R — » oo ) and very small at x = xj (i. e. 
R = R min)> with some intermediate value R| , at x = x m . This assumes that 
the interface in the evaporator is in the position shown in Figure 5 (a). 

Integrating Equation (3-6) between R = R ^ (x = 0) and R = Rj (x = x m ) yields 


2 <r 
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(3-12) 


Integrating Equation (3-11) between R = Rj (x = x m ) and R = R m i n ( x = x^) 
yields 
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(3-13) 
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Addition of Equations (3-12) and (3-13) yields 


2 <r 
Rmin 


K P T a — 

1 Lh vi* p : 


X Xrp 

m 1 


S-x-j- P T sin e = 0 

or L 


(3-14) 


There are no momentum terms in Equation (3-14) which is expected since the 
liquid velocity is zero at x = 0 and x = xt . 


As shown in Equation (47) 


R . = 
min 


V g *m 


) 


WR 


Substituting R . into Equation (3-14) and rearranging terms yields an expression 
for the maximum total heat input that can be transferred to a vapor-chamber fin 
as a function of geometry, wick properties and fluid properties. 



(3-15) 


Neglecting the evaporator section (i. e x T = x^) and solving for the maximum total 
heat transferred to a fin in a horizontal position Equation (3-15) reduces to 
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max 



L h VL 


P- 


<r 



which is identical to Equation (49) in Section in D. 


(3-16) 
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